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Reference  is  made  to  the  mid-term  report.  In  addition,  the  simulation 
study  I has  been  completed  by  a more  general  model  in  simulation  study  II. 

A detailed  report  of  115  pages  is  added  to  this  summary. 

The  new  analysis  of  geodetic  applications  of  inertial  navigation  systems  started 
from  the  basic  equation 

» 

"observed  acceleration  “ 

Inertial  acceleration  minus  gravitation" 

which  can  be  solved  for  inertial  acceleration  if  we  know  gravitation,  e.g.  from 
integrated  gravity  gradient  measurements.  The  naive  approach  is  as  following: 
Measure  the  coordinates  of  the  apparent  acceleration  vector,  e.g.  Ax,  Ay,  Az, 
and  the  coordinates  of  the  tensor  of  gravity  gradients,  e.g.  Gxx.  Gxy,  Gxz» 

Gyy,  GyZ,  Gzz . Then  compute  at  point  one 

Axi  + GXo  + Gxx0  (Xl  “ *o)  + Gxyo  (Y1  “ Yo)  + Gxz0  (Zl  “ Z0)  “ 

At“2(X2  - 2X1  + XD) 

Ay j , AZj  analogous 

where  At  indicates  the  time  interval  of  measurements;  the  right  side  approxi- 
mates Inertial  acceleration  by  Stirling’s  formula. 

What  we  have  done  at  point  one  can  be  done  at  any  point  such  that 
*i+l  " ^[Axa  -Wxi  + GxXi.iOCi  - + GJCy1_1  (Yi  - Yt_p 

+ GXzj_i(zi  - Zl.p]  + quantization  errors 
Yi+1>  zi+l  analogous 
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are  the  unknown  coordinates  of  points  with  respect  to  inertial  space  we 
like  to  know.  The  set-up  is  in  terms  of  an  initial  value  problem  since 
we  have  to  know  (X0,  Y0,  ZQ)  (Xj,  Yj,  Zj},  the  coordinates  of  the  initial 
position  vectors,  and  (G^,  GyQ,  GZqJ,  the  coordinates  of  the  initial 
gravity  vector.  In  addition,  the  Laplacean 

Gxx  + Gyy  + Gzz  * 0 

is  a physical  condition  if  we  measure  outside  the  masses  (surfaces  in- 
cluded) . Finally,  a transformation  of  inertial  coordinates  into  terres- 
trial ones  has  to  be  performed. 

Simulation  I is  an  error  budget  study  by  19  parameters  of  an  inertial 


system: 

(i) 

time  interval 

(ID 

Initial  positions 

(111) 

initial  gravity 

(iv) 

varying  acceleration 

(v) 

varying  gravity  gradients 

The  input-output  results  are  given  on  pages  46  - 56  of  the  detailed  report. 
Pages  61  - 67  are  a study  of  the  influence  of  quantization  errors.  This  is 
a summary: 

The  gradiometer  accuracy  of  existing  systems  is  by  far  sufficient  for  inertial 
navigation  applications;  the  main  error  budget  is  due  to  the  accelerometer 
accuracy . 

Simulation  II  is  an  error  budget  study  by  36  parameters  of  an  inertial 


system: 

(1) 

time  interval 

(ID 

initial  positions 

(ill) 

initial  gravity 

(iv) 

varying  acceleration 

(v) 

varying  gravity  gradients 

(vi) 

accelerometer  bias 

(vii) 

accelerometer  random  uncertainty 

(viii) 

accelerometer  non-orthogonality 

(lx) 

initial  misalignment  angles 

(x) 

accelerometer  scale  factor  uncertainty 

The  input-output  results  are  given  on  pages  96  - 108  of  the  detailed  report. 


March  26,  1979 
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Abstract 


A gradiameter-aided  inertial  navigation  system  is  theoretically  and  statistically  analysed  tc 
estimate  its  abilities  to  monitore  geocentric  cartesian  coordinates.  Having  discussed  the  inertial 
instrumental  units  used  on  the  moving  platform  and  several  reference  coordinate  frames  applicable 
in  all  navigation  systems,  studies  on  the  severe  problem  of  the  separability  of  the  gravity 
gradients  from  the  inertial  disturbances  are  carried  out.  Simulation  1 presents  how  well  the 
aided  navigation  system  can  produce  inertial  coordinates  and  how  the  newcomers  of  the  inertial  in- 
strumentation, the  gradiometers,  perform  on-board  the  moving  vehicle.  Quar.ti ration  error  studies 
are  also  analysed  and  presented  for  such  a system.  Simulation  II  includes  besides  tne  detailed 
analysis  of  the  accelerometer  and  gradicmeter  error  models  used,  the  abilities  of  the  system  to 
estimate  geocentric  coordinates.  Multipoint  statistical  analysis  for  the  approximated  inertial 
acceleration  components  shows  that  the  navigation  system  under  consideration  behaves  setter  as 
closer  the  reality  is  approached. 


lusanmenfassur.g 

Lin  cradiometer-unterstutrtes  inertiales  Navigations system  w. rd  theoretisch  und  statistiscl  analy- 
siert,  or.  seine  Moglichkeiten  abzuschStzen,  geozentrische  cartesische  Koordinaten  zu  ermitteln. 
Nach  einer  Disfcussion  der  ir.ertialen  Instrunenten-Einheiten  aui  der  bewegten  Plattfonr  und  eini- 
gtr  Bezugssysteme,  die  in  alien  Navigationssystenen  angewandt  werden,  wird  das  Problem  der  Trenn- 
barkeit  der  Schweregradienten  von  den  Inertialstbrungen  untersucht.  Simulation  1 zeigt,  mit  Wel- 
ches Giitc  das  untcrstiitztc  Navigationssystem  inertiale  Kocrdincten  lie-fern  kann  und  wie  die  Neu- 
linge  unter  den  Inertialgeraten,  die  Gradiometer,  sich  an  Bord  des  bewegter.  Farrzeuges  verbal  ten. 
Auch  der  Einflufl  von  Qiantisierungsfehier.i  wird  untersucht  und  fUr  soldi  em  Syster..  present iert. 
Simulation  II  enth&lt  neben  einer  detaillierten  Analyse  der  benutzten  Fehiermodelle  fur  Accelero- 
meter und  Gradiometer  auch  die  Moglichkeiten  des  Systems  zur  Schatzung  geozer.tr ischer  Kcirdinaten. 
Fine  statistische  Ntehrpunkt-Analyse  der  nAhe rungs we i se  inertialer.  Deschieunigungskcmponen'en 
zeigt,  da/1  das  betrachtetc  Navigationssystem  sich  unso  bessei  verhfllt,  jc  nAher  man  der  Real.’tiit 
taint. 
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0.  Introduction 


Paris,  May  1914.  Lawrence  Sperry  is  one  of  the  entrants  for  a prize  of  fifty 
thousand  francs,  the  object  of  the  competition  being  what  one  judge  sees  and 
describes:"the  airplane  is  already  in  flight  and  the  mechanic  rises,  leaves 
his  seat  and  without  fear  goes  from  the  cabin  to  the  wings  and  returns  back. 

At  the  same  moment  the  pilot,  Lawrence  Sperry,  lifts  his  arms  and  leaves  the 
airplane  to  continue  its  flight  without  guidance  with  speed  of  about  90km/n". 

He  finally  won  and  the  first  stabilized  platform  was  already  introduced  for 
airplane  stabilization. 

Navigation  is  perhaps  connected  with  the  first  human  activities  or.  the  earth 
because  it  is  the  "art ' of  obtaining  with  measurements  velocity  and  position 
of  a moving  object.  Guidance  has  tc  be  always  distinguished  from  navigation, 
for  it  is  the  process  of  generating  motion  correction  commands  to  a rxvitg 
object  such  as  it  succeeds  in  its  mission,  it  is  understood  tnat  guidance 
includes  navigation  but  not  vice-versa. 

Guidance  techniques  could  be  considered  as  an  extension  of  human  being’s  natu- 
ral senses  of  seeing,  smelling,  hearing,  feeling,  memory  and  deduction.  In  ca- 
ses where  the  guidance  problem  is  not  inside  the  man’s  abilities,  then  a "de- 
vice" could  upgrade  one  or  more  of  his  natural  senses.  For  example,  man’s  mea- 
suring abilities  could  be  augmented  by  using  electronic  measuring  devices  or 
computers  which  can  increase  his  accuracy  and  speed  of  deduction.  Going  a lit- 
tle further  where  the  guidance  problem  is  very  complex  or  the  human  presence 
is  impossible  (e.g.  mission  leading  to  destructive  ten.inatior) , then  automatic 
guidance  is  to  be  introduced. 

Newton’s  law  of  mechanical  inertia  is  by  far  the  basic  law  governing  all  desi- 
red properties  of  guidance  systems.  Under  this  law  all  particles  with  mass  will 
exert  reaction  to  the  applied  acceleration,  which  it  equal  in  magnitude  and  op- 
posite in  direction,  the  reaction  not  being  depenucr.t  on  contacts  with  the  envi 
ronment . Constructing  self-contained  instruments  dependent  upon  inertial  effects 
it  is  therefore  natural  to  refer  to  them  as  inertial  guidance  systems  and  in  ca 
se  of  navigation  as  inertial  navigation  systems. 

An  inertial  navigation  system  generally  contains  four  basic  elements':a)an  accele 
rometer  b)ai  attitude  ref srencec) a computer  and  a)a  clock.  An  accelerometer  is  a 
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device  which  measures  the  non-gravitational  acceleration  experienced  upon 
its  case.  But,  since  the  case  is  hard -mounted  on  the  moving  vehicle,  the 
measured  acceleration  is  also  the  acceleration  of  the  vehicle.  It  is  under- 
stood that  since  a vector  quantity  has  three  components , then  three  accelero- 
meters', orthogonally  mounted,  measure  the  non-gravitational  acceleration  re- 
solved in  their  sensitive  or  input  axes.  The  orientation  of  these  axes  is  in 
a manner  best-suited  to  meet  the  system's  requirements. 

The  attitude  reference  is  that  part  of  the  inertial  navigation  instrumenta- 
tion which  either  stabilizes  or  commands  the  platform  frame  relative  to  an 
inertial  or  rotating  frame  respectively.  Gyroscopes  are  always  called  to 
instrument  the  attitude  reference  and  considered  to  be  the  most  indispensable 
and  critical  unit  on-board. 

The  conputer  solves  the  fundamental  equation  of  inertial  navigation  to  give 
the  velocity  and  position  estimates  of  the  vehicle  being  navigated.  The  ma- 
thematical procedure  to  be  followed  depends  on  the  system  requirements  as  veil  as 
the  on-board  instrumentation.  If,  for  example,  the  accelerometers  used  are 
of  the  integrated  type,  then  the  conputer  does  not  perform  integrations  as 
far  as  the  accelerometer  signal  is  concerned. 

The  clock  generally  gives  the  time  instant  of  the  measurements  performed  on- 
board. In  cases  the  navigation  takes  place  relative  to  an  earth  rotating  fra- 
me, then  the  clock  establishes  the  location  or  orientation  of  that  frame. 

An  inertial  navigation  system  may  be  abstractly  considered  as  a black  box. 

The  input  to  this  box  is  apparent  acceleration  which  contains  relative  ve- 
hicle acceleration,  gravitational  acceleration  and  accelerations  due  to  the 
rotation  of  the  frame  of  the  black  box  relative  to  the  inertial  space.  In- 
side the  black  box  manipulations  are  carried  out  and  the  output  is  finally 
the  instantaneous  velocity  and  position  of  the  moving  vehicle  or  better  of 
the  black  box.  Needless  to  say,  that  the  manipulations  contain  errors  and 
therefore  the  output  is  incorrectly  indicated.  The  differences  between  the 
actual  velocity  and  position  and  their  output  counterparts  are  then  stati- 
stically analysed  to  yield  the  error  budget  of  the  navigation  system (black 
box) . 

It  could  be  supposed  that  inertial  navigation  is  acconplished  by  measuring 
the  apparent  acceleration  of  the  moving  vehicle  and  then  integrating  it 
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doubly  in  time  to  estimate  position.  But  due  to  Einstein’s  principle  of 
equivalence,  the  gravitational  acceleration  and  the  non-gravitational  one 
are  manifestations  of  the  same  basic  physical  phenomenon.  Consequently, 
gravitational  information  is  certainly  needed  on-board  in  order  to  solve 
the  navigation  equation  and  there  is  no  way  out  of  it. 

Tracing  the  historical  developments  of  inertial  navigation,  we  see  how  the 
scientists  treated  the  above  postulation  due  to  the  lack  of  actual  gravita- 
tional information.  A reference  field  was  selected,  the  gravitational  acce- 
leration was  confuted (grossly  approximated) at  the  instantaneous  vehicle  lo- 
cation and  then  substructed  from  the  apparent  acceleration  measurements  to 
go  to  the  solution  of  the  final  equation  used.  Unfortunately,  the  same  pro- 
cedure is  still  in  use  nowadays. 

As  early  as  in  '950,  Lundberg  constructed  and  tested  the  first  gravity  gra- 
diometer.  The  instrument  was  composed  of  two  vertically  suspended  masses  and 
had  the  ability  to  sense  the  sign  of  the  first  vertical  derivative  of  gravi- 
ty. The  gradiometer  was  heavily  tested  in  North  America, Europe  and  West  Atti- 
ca, but  nowadays  is  completely  forgotten  for  not  known  reasons. 

During  1959-62,  the  Lockwood  Company  in  Toronto  introduced  its  first  prototy- 
pe gradiometer  sensing  the  vertical  gradient  of  gravity  with  an  accuracy  of 
_q  _? 

100Ebtvos(1E=10  'sec  “) , but  the  device  was  very  sensitive  to  the  airborne 
dynamic  environment. 

After  these  first  attempts  came  the  era  of  the  fi^st  generation  of  gravity 
gradio meters.  The  Hughes,  the  M.I.T,  the  Bell  Aerospace  gradiometers  are  a 
sample  of  very  refined  and  tested  gradiometers  designed  to  be  used  in  air- 
borne gradiometry.  The  feasibility  studies  have  proved  so  far  tr?t  an  accu- 
racy of  IE  or  better  is  to  be  expected  in  the  very  near  future  and  especial- 
ly for  the  Hughes  gradiometer  for  which  five  days  of  continuous  gathering 
data  of  the  earth’s  gravity  field  would  be  enough  to  map  it  completely. 

Savet  in  one  of  his  papers  writes: "altogether,  it  appears  that  there  is  no  ^ 
clear-cut  advantage  in  using  an  existing  or  feasible  gradiometer  or,  for 
that  matter,  a pair  of  accelerometers" (Savet,  1970).  Eight  years  later,  we 
read  something  breathtaking: "in  the  University  town  of  Nancy  in  France 
George  Delamare  declares :« the  key  to  my  operation  is  a tiny  electrode  inp- 
lanted  in  each  of  the  leg  muscles  which  transmits  a computerized  electro- 
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stimulation  in  response  to  a peripheral  data  system  based  upon  micro  acce- 
lerometers and  inclinometers  to  appreciate  the  patient (paraplegic) space 
mission.  Experiments  carried  out  in  the  NASA  space  programme  could  be  ap- 
plied to  the  problem  of  controlling  the  balance  of  the  patient  when  mo- 
bile>>"(Time  magazine,  March  1978). 

It  seems  that  sometimes  non-geodesists  appreciate  the  geodetic  tools  better. 

In  view  of  the  era  which  the  gradiometers  promise  and  the  fact  that  there 
is  rot  such  a gross  exaggeration  as  to  assume  a spherical  earth  gravity 
field(!),  we  motivate  the  present  analysis  to  put  on-board  the  moving  ve- 
hicle a number  of  gradiometers  to  measure  the  earth’s  gravity  field.  One 
could  immediately  assert  that  the  accuracy  of  the  inertial  navigation  sy- 
stems is  slightly  better  than  one  mile/hour  flight  and  it  might  be  consi- 
dered as  satisfactory.  But,  again,  we  do  not  lay  the  problem  on  the  desi- 
red or  obtainable  accuracy.  Perhaps  we  do  not  think  so  much  of  the  opera- 
tional point  of  view,  at  least,  in  the  very  beginning.  The  navigation  sy- 
stems need,  as  we  believe,  a theoretical  injection  which  inevitably  comes 
from  the  gradiometer  implementation. 

Consequently,  we  start  the  whole  analysis  fresh  from  the  beginning.  We  avoid 
one  of  the  current  techniques,  that  is,  first  make  the  assumptions (so  bias 
the  system)  and  then  obtain  what  it  might  be  expected  .-good  results  like  the 
assumptions.  We  go  the  other  way  round.  Lay  the  fundaments  rigorously  and 
the  time  for  the  assumptions  will  come  in  order  to  present  some  indicated 
numbers  of  how  well  the  system  behaves.  Perhaps,  there  is  a reconciliation 
between  the  two  approaches. 

The  newcomers  in  the  inertial  navigation  instrumentation  technology,  the  gra 
diometers,  promise  many  applications.  As  we  all  know,  one  of  the  prime  goals 
in  geodesy  is  the  determination  of  the  earth’s  gravity  field.  This  can  be 
straightforwardly  acconplished  by  employing  a moving  gravity  gradiometer 
which  could  furnish  even  in  a very  short  period  the  desired  gravity  data. 
Other  applications  of  a gradiometer-aided  inertial  navigation  system  are: 

a) vertical  indication 

b)  vertical  deflection  indication 

c) geoid  height  indication 

just  to  mention  only  a few  of  them.  It  is  therefore  concluded  that  gravity 
gradiometers  offer  valuable  applications  in  geodetic  science. 

I 
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The  inertial  navigation  system  we  analyse  is  a gradiometer-aided  system. 
Whatever  comes  as  an  input  into  the  black  box  is  actually  measured  on-board . 

Certain  mathematical  manipulation  takes  place  inside  the  box  and  the  output, 
what  we  are  particularly  interested  in,  is  instantaneous  vehicle  coordinates 
with  respect  to  a frame  rotating  with  the  earth. 

Section  1 deals  with  the  inertial  navigation  instrumental  units.  Our  plat- 
form is  enriched  by  three  spherical  gradiometers  developed  and  tested  in 
the  M.I.T.  Three  gyros  command  the  platform  to  rotate  with  the  earth’s  rate 
and  the  system  has  the  capability  to  actually  measure  and  then  feed  the  he 
on-board  conputer  with  all  what  it  needs : apparent  acceleration  components 
and  the  whole  gravity  gradient  tensor. 

In  section  2,  the  fundamental  equation  of  inertial  navigation  for  an  earth- 
bound  region  is  derived.  Since  we  are  only  interested  in  examining  the  capa- 
bilities of  a gradiometer-aided  inertial  navigation  system  for  terrestrial 
navigation,  the  gravitational  fields  of  all  attracting  masses  but  the  earth’s 
are  excluded. 

The  description  of  all  common  coordinate  frames  used  in  inertial  navigation 
systems  is  presented  in  section  3.  As  scon  as  the  notion  of  these  frames  is 
completely  understood,  then  one  could  really  have  the  flexibility  to  intro- 
duce the  most  general  distortions  a coordinate  frame  can  undergo.  Since  our 
navigation  system  estimates  geocentric  coordinates  of  the  mass  centre  of 
the  platform  and  it  is  impossible  to  "curve  in"  all  on-board  frames  on  tnat 
point,  due  to  the  actual  dimensions  of  the  units,  it  is  assumed  that  the  re- 
sulted centrifugal  accelerations (each  for  each  instrumental  unit  mass  centre) 
are  negligible  small  quantities. 

Thinking  in  terms  of  Einstein’s  principle  of  equivalence,  it  seems  to  be  ho- 
peless to  separate  gravitational  from  inertial  effects.  Studies  on  this  se- 
parability are  carried  out  in  section  4.  It  is  concluded  that  gradiomcter; 
can  measure  in  a dynamic  environment  only  gravity  gradients  if  and  only  if 
they  are  inertially  stabilised.  Since  our  platform  is  inertially  rotating, 
a second  platform  is  introduced  to  accommodate  only  the  gradiometer  measu- 
rement unit  stabilized  with  respect  to  an  inertial  frame. 

An  extensive  literature  exists  on  simulation  studies  of  inertial  navigation 
systems,  Depending  on  the  system  used  and  the  assumptions  set,  different  re- 
sults have  been  drawn  so  far.  The  performance  of  the  inertial  instrumentation 
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is  considered, nowadays,  to  be  satisfactory,  but  the  systems  still  gather  an 
appreciable  error  budget.  Updating  the  system  from  time  to  time  , the  errors 
are  reduced  but  they  are  too  far  from  being  eliminated.  Since  there  is  not 
any  actual  in  flight  gradiometer,  it  is  assumed  and  sometimes  believed  that 
the  approximation  of  the  earth's  gravity  field  has  the  worst  contribution 
into  the  system’s  accuracy.  For  that  reason,  we  try  to  identify  with  a 
simple  example  the  function  which  drives  the  errors  in  an  inertial  naviga- 
tion system.  Section  5 contains  a simulation  to  this  direction  not  inclu- 
ding error  models  for  the  instrumental  units.  Finally,  a space  traverse  is 
computed  to  see  the  behaviour  of  such  an  aided  system. 

Section  6 examines  the  quantization  of  errors  for  the  accelerometers  and 
gradiometers.  Assuming  that  each  time  the  instruments  are  read  a quantiza- 
tion error  is  present,  a general  formula  is  given  to  estimate  the  position 
error  due  to  the  instrumental  truncations  treated  as  stochastic  quantities. 

A general  error  model  is  analysed  in  section  7.  We  go  to  the  far  end  of 
each  measurement  unit  taking  into  account  what  they  actually  measure  and 
considering  the  possible  types  on  instrumental  misalignment  and  non-ortho- 
gonality. The  errors  which  contribute  more  than  90 % of  the  whole  errcr  budget 
are  paid  attention  and  included  in  the  general  model.  On  the  final  navigation 
equation  suitably  approximated,  the  second  simulation  studies  are  performed 
to  test  the  capabilities  of  a gradiometer-aided  inertial  navigation  system. 


I 
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1 . Inertial  measurement  units 

1.1  Accelerometers 

1.1.1  General  considerations 

Just  a glance  at  the  fundamental  equation  of  inertial  navi gation(eq. (2.1 1)) 
shows  how  critical  is  the  acceleration  contribution  to  the  navigation  pro- 
blem. As  a matter  of  fact,  inertial  navigation  can  be  principally  accom- 
plished by  measuring  only  the  apparent  acceleration  on-board  a moving  ve- 
hicle. The  gravity  field  compensation  could  be  computed  by  employing  a re- 
ference field,  such  as  an  ellipsoidal  one,  and  using  approximate  position 
values,  so  as  the  gravity  components  at  the  instrumentation  location  are 
then  to  be  estimated. 

The  lumber  of  the  on-board  accelerometers  depends  on  the  particular  pro- 
blem considered.  Generally  speaking,  when  we  navigate  in  the  three  dimen- 
sional space,  three  accelerometers  are  used  in  order  to  sense  all  three 
apparent  acceleration  components.  In  case  of  a cruise  aircraft,  the  verti- 
cal accelerometer  could  be  substituted  by  a barometer  or  an  altimeter  re- 
ducing the  number  to  two. 

The  accelerometers  provide  their  measurements  in  a frame  which  they  are 
constrained  to  follow  and  it  is  furnished  by  the  gyros.  In  our  analysis, 
three  accelerometers  are  set  at  the  platform’s  mass  centre  such  as  to  con- 
struct an  orthogonal  frame  and  measure  the  apparent  acceleration  vector 
resolved  in  their  rotating (with  earth’s  rate)  axes. 

1.1.2  Operational,  principles 

Let  us  now  see  how  an  ideal  accelerometer  operates  and  in  order  to  moti- 
vate our  discussion,  consider  an  ideal  spring-mass  accelerometer  shown  in 
Fig.(t).  The  instrument  consists  of: 

a)  the  case 

b)  the  proof  mass  and 

c)  the  damping  spring 

As  the  accelerometer  follows  the  motion  of  the  vehicle,  the  apparent  acce- 
leration acts  on  the  proof  mass  as  well  as  the  case.  The  proof  mass  extends 
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^ Input  axis 


Fig. (1 ) : An  ideal  spring-mass  accelerometer 

the  spring  and  the  displacement  x is  measured  against  a properly  divided 
scale.  Generally,  the  equation  governing  the  spring-mass  accelerometer 
operation  reads 


(i.i.i) 


mx*d*+kx 


where  A the  apparent  acceleration,  x the  scale  reading,  m the  proof  mass, 
d the  damping  factor  and  k the  spring  constant.  Having  got  the  x scale  rea- 
ding and  known  the  mentioned  constants,  the  apparent  acceleration  can  be 
then  computed  using  eq. (1.1.1).  Needless  to  say,  that  the  apparent  accele- 
ration is  sensed  along  the  input  accelerometer  axis  and  therefore  three 
such  accelerometers  orthogonally  mounted  can  fully  estimate  the  apparent 
acceleration  vector. 

Finally,  we  remark  that  an  ideal  accelerometer  is  nothing  else  but  an  ideal 
gravimeter  and  since  a moving  gravimeter  cannot  discriminate  between  gravi- 
tational and  non-gravitational  forces,  the  accelerometer  output  is  a mix- 
ture of  these  two  force  fields  and  as  such  it  must  be  gravimetrically  com- 
pensated . 

1.1.3  Error  model 


As  we  said  before,  three  accelerometers  can  measure  the  apparent  acceleration 
vector  and  their  output  contributes  to  the  estimation  of  the  inertial ly  re- 
ferenced acceleration.  But,  since  it  is  inpossible  to  construct  an  orthogonal 
accelerometer  frame,  the  apparent  acceleration  is  sensed  along  a quasi-ortho- 
gonal  frame  and  care  should  be  taken  to  compensate  foi  accelerometer  non- 
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orthogonality.  The  transformation  of  the  apparent  acceleration  vector 
from  the  quasi-orthogonal  to  the  ideal  orthogonal  accelerometer  frame 
is  a small  angle  transformation  and  it  is  considered  in  full  detail  in 
section  7 . 1 . Furthermore , each  accelerometer  has  its  own  bias  and  scale 
factor  uncertainty  which  falsify  the  sensed  acceleration.  Taking  all 
these  errors  into  account,  the  accelerometer  signal  is  corrected  and 
refers  to  the  ideal  accelerometer  frame  ready  to  be  transformed  to  any 
desired  coordinate  frame  in  which  we  like  to  solve  the  fundamental  equ- 
ation of  inertial  navigation. 

1 . 2 Gradiometers 
i . 2 . i General  considerations 

Right  frcr.  the  beginning  of  the  application  of  the  inertial  navigation 
systems,  it  was  fully  understooe  that  a serious  contributing  factor  in 
their  error  budget  was  the  gravity  field  compensation  in  the  accelero- 
meters measurements.  At  that  time,  there  was  no  space  for  envisioning 
an  instrument  mounted  on-boara  a moving  vehicle  and  having  the  capabi- 
lity of  directly  measuring  the  gravity  field.  As  the  years  went  by  a ri 
the  inertial  instrumentation  and  guidance  systems  reached  such  a tre- 
mendous qualification  and  performance,  it  was  quite  evident  that  a re- 
search for  the  first  generation  of  moving  gravity  instrumentation  was 
inevitable.  Sc  far,  the  earth’s  gravity  field  was  approximated  by  the 
gravity  field  of  an  ellipsoid  of  revolution.  Giver,  approximate  position 
values,  the  gravity  field  is  computed  at  the  vehicle’s  instantaneous 
location  using  well-known  closed  formulas.  Having  granted  by  thr  gra- 
dioneter  instrumentation,  we  feel  a little  disturbed  to  continue  appro- 
ximating the  reality  as  far  as  the  gravity  field  is  concerned.  Of  course, 
the  motivation  for  studying  and  constructing  gradiometers  was  not, at 
least  primarily,  to  aid  inertial  navigation  systems.  But  since  gravity 
is  a critical  inertial  parameter,  we  think  that  it  is  really  a good >q 
chance  to  have  on-board  a gradiometer,  or  a number  of  them,  in  an  effort 
to  actually  measure  the  gravity  field  in  which  the  navigation  takes  pla- 
ce. 

Nowadays,  the  gradiometers  are  still  in  the  feasibility  study  phase  except 
some  on-board  tests (see  P.Hood  and  H. Ward, 1969, p. 98) , but  their  performance 
promises  very  quick  on-board  inplementation(Trageser,  1975) . Consequently, 
it  is  reasonable  to  assime  that  under  the  assertions  of  the  mechanical 
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engineers,  gradiometer  applicability  is  warranted  in  a serious  effort  to 
aid  autonomous  inertial  navigation  systems. 

1.2.2  Operational  principles 

A gradiometer  measures  the  changes  of  the  gravity  conponents  with  respect  to 
the  displacement, that  is, gravity  gradients.  The  advantages  of  a gradiometer 
arereally  a lot.  Moving  base  capability,  no  E6tv5s  correction,  no  terrain 
correction,  just  to  mention  a few  of  them.  Among  the  most  interesting  ones 
are: 

a)  the  Bell  Aerospace  gradiometer 

b) the  Hughes  gradiometer  and 

c) the  M.I.T.  gradiometer 

For  a description  see (Williams, 1975) . In  1915,  the  intelligencv  of  Rolan 
von  Eotvbs  created  the  torsion  balance  bearing  his  name  and  measuring  cer- 
tain gravity  gradients.  Without  exception,  its  principles  apply  to  the  up- 
to-date  gradiometers . Consider  two  equal  proof  masses  connected  with  an 
axle  and  supported  at  a flexure  point  between  them.  When  that  primitive 
instrument  passes  over  a mass  anomaly,  different (in  the  differential  sense) 
attraction  forces  are  exerted  upon  the  proof  masses  due  to  the  difference 
in  the  distance  between  them  and  the  mass  anomaly.  Consequently,  the  instru- 
ment changes  its  position  about  the  flexure  point  and  if  this  change (ro- 
tation) is  appropriately  sensed,  then  gravity  gradients  could  be  measured. 

Of  course,  the  measurement  of  a gradient  is  very  complicated  due  to  the 
advanced  electronics  involved. 

In  what  follows,  we  shall  try  to  explain  briefly  the  spherical  gradiometer 
developed  in  the  Massachusetts  Institute  of  Technology (M. I .T.)  which  is 
used  as  our  on-board  gradiometer.  Since  all  the  known  gradiometers  are'uhdeT 
laboratory  tests .there  is  not  any  immediate  justification  why  we  prefer  that 
instrument  to  the  others.  But  we  believe,  at  least  from  the  given  accuracies 
(Trageser,1975,  Forward,  1974,  Williams,  1975),  that  the  M.I.T.  spherical 
gradiometer  possesses  certain  advantages  in  comparison  with  the  others  e.g. 
structural  stability,  immunity  to  platform  jitter  rectification  effects  etc. 

Let  us  therefore  study  the  operational  principles  of  the  spherical  gradiome- 
ter. The  instrument  consists  of  (see  Fig. (2)): 

a)  the  float  and 

b)  the  housing  of  the  float. 
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Fig.  (2):  Spherical  gradiometer  float 

Between  the  float  and  the  housing,  a small  gap  exists  which  is  filled  by 
a special  fluid  supporting  the  float  inside  the  stabilized  housing.  When 
the  instrument  passes  over  a mass  anomaly,  the  float  is  constrained  tc  ro- 
tate inside  the  housing.  The  moment  about  a defined  axis  is  measured  by 
the  electronics  of  the  housing  and  is  applied  back  as  a restoration  torque 
in  order  to  bring  the  float  to  its  original  position.  For  example,  the  mo- 
ment measured  about  the  X axis  (Fig.  (2))  gives  the  gy-,  gravity  gradient 
and  the  moment  about  the  Y axis  the  g^.~  one.  Therefore,  each  instrument 
measures  two  gradients (another  great  advantage)  and  three  of  them  mounted 
on  a platform  can  provide  us  the  whole  gravity  gradient  tensor. 

Feasibility  studies  on  the  spherical  gradiometer  support  the  expectation 
that  IE  or  better  could  be  attained  in  the  near  future (Trageser,  ’975) 

1.2.3  Error  model 

In  trying  to  construct  a general  gradiometer  error  model,  three  reference  <■’* 
frames  are  employed  in  conduction  with  the  gradiometer  float: 

a)  the  electronic  frame,  which  senses  the  small  float  rotations (torques) 

b)  the  float  or  misaligned  gradiometer  measurement  unit  frame,  vchich 

is  fixed  with  respect  to  the  float  principal  moments  of  inertia  axes  and 
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c)  the  ideal  gradiometer  measurement  unit  frame  to  which  the  gravity  gra- 
dients refer. 

In  a general  error  model,  the  first  two  reference  frames  should  be  consi- 
dered as  a quasi -orthogonal  ones  and  the  gravity  gradients  signal  must  be 
accordingly  corrected  using  the  involved  misalignment  angles.  Furthermore, 
gradiometers  are  instruments  which  inertial ly  measure  gravity  gradients. 
Since  our  gradiometer  measurement  unit  frame  is  designed  to  have  the  same 
orientation  as  an  earth  centered  quasi-inertial  frame,  a general  misalign- 
ment of  the  former  frame  with  respect  to  the  latter  one  should  be  taken 
into  account.  When  all  these  misalignments  induced  gravity  gradients  errors 
are  considered,  then  the  gravity  gradients  signal  can  be  transformed  to 
the  desired  coordinate  frame  in  which  we  like  to  solve  the  fundamental  equa- 
tion of  inertial  navigation. 

1 .3.  Gyroscopes 

1.3.1.  General  considerations 

It  is  well-known  that  inertial  navigation  reaches  its  aim  in  determining 
velocity  and  position  with  respect  to  a reference  frame  by  the  implementa- 
tion of  three  mutually  perpendicular  mounted  accelerometers  which  measure 
the  components  of  the  apparent  acceleration  vector  resolved  in  their  input 
or  sensitive  axes.  It  is  noting  by  passing,  that  these  measurements  must  be 
gravitationally  compensated.  Let  us  now  suppose  that  we  like  to  measure  the 
apparent  acceleration  in  a moving  vehicle  in  order  to  determine,  say,  posi- 
tion. For  that  purpose,  three  accelerometers  are  mounted  on-board  the  mo- 
ving vehicle  having  their  sensitive  axes  to  a known  orientation  relative 
to  a reference  frame  with  respect  to  which  we  navigate.  As  soon  as  the  ve- 
hicle starts  moving  and  begins  pitching,  rolling  and  yawing,  then  it  is 
absolutely  inpossible  for  the  accelerometer  axes  to  preserve  their  original 
orientation.  Consequently,  their  output  cannot  be  used  as  an  input  for  sol- 
ving the  fundamental  equation  of  inertial  navigation.  A special  device  must 
be  inplemented  in  order  to  ’’dump  out”  the  time-iik.2  accelerometer  frame  miso 
rientations.  This  device  should  have  the  capability  of  either  to  command  the 
accelerometer  frame (or  for  that  matter  the  platform  frame)  to  rotate, in  case 
the  navigation  frame  rotates  or  to  stabilize  it  to  a desired  orientation, 
in  case  of  an  inertially  non-rotating  reference  frame.  This  can  be  accomp.- 
lished,  as  we  will  see,  by  employing  three  single-degree-of  freedom  gyros 
or  two  twp-degree-of  freedom  which  can  cover  the  three  possible  degrees  of 
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freedom  in  space. 

In  what  follows,  we  shall  try  to  give  briefly  the  gyro  functional  chara- 
cteristics and  their  error  model  which  is  included  in  the  simulation  stu- 
dies. 

1.3.2  Operational  principles 

Let  us  suppose  that  a solid  body  rotates  with  angular  velocity  u along 
an  axis  of  symmetry-  (Fig.  (3a)).  Newton’s  second  law  in  rotational  form 
states  that  in  an  inertial  frame  of  reference: 


where  f represents  the  applies  torque,  H the  angular  momentum  and  I the 
moment  of  inertia  of  the  rotating  solid  mass.  Now,  let  mount  this  simple 
device  on  a case  with  a single  axis  gimlal  as  in  (Fig. (3b)).  This  over- 
simplified mechanization  depicts  the  principle  of  the  single-degree-of- 
freedom  gyro (SLF)  owing  its  name  co  the  single  gimbal  suspension.  A SI3F 
gyro  is  composed  of: 


Fig. (3):  Single-degree-of -freedom  gyro. 


a)  the  gyro  case  or  housing 

b)  the  gyro  float  and 

c)  the  gyro  rotor 

From  eq. (1.3.1),  it  is  understood  that  in  absence  of  applied  torques,  its 
integration  gives 

ft  * constant 

and  this  manifests  the  most  fundamental  property  of  a gyro:when  no  external 
torque  is  present,  then  the  direction  of  the  angular  momentum  vector  pre- 
serves its  orientation  in  space. 


In  order  to  practically  utilize  the  property  of  the  constancy  of  the  angu- 
lar momentum  orientation,  the  H vector  ought  to  be  in  the  direction  of  one 
axis  of  the  rotating  body  symmetry.  If  it  is  not  the  case,  then  the  angular 
momentum  could  be  resolved  in  three  components  as: 


(1.3.2) 
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where  u>s  the  rotation  along  the  axis  of  symmetry  and  ,u>2  the  rotations 
along  two (any)  perpendicular  axes.  Denoting  the  angle  between  ft  and  the 
axis  of  symmetry  by  we  note  that  this  is  a constant  angle,  since  it 
depends  on  the  I’s  and  w’s  which  are  constant  quantities  for  a solid  body 
rotating  with  constant  speed.  But  due  to  the  principle  that  the  angular 
momentum  vector  preserves  its  orientation  in  space,  then  it  is  concluded 
that  the  axis  of  symmetry  must  move  otherwise  it  would  coincide  with  ft. 
Consequently,  the  axis  of  symmetry  generates  a cone  around  H with  apex 
angle  2\i>  and  we  refer  to  it  as  the  (free)  gyroscopic  nutation  phenomenon. 


Suppose  now  that  a torque  is  applied  perpendicular  to  the  angular  momentum 
vector  H.  In  this  case,  H rotates  with  angular  velocity  3 transverse  to 
both  torque  and  ft  vectors  and  it  is  such  that  the  applied  torque  is  equal 
to  the  cross  product  of  ft  and  the  rotation  vector  1§.  This  rotation  of  the 
angular  momentum  vector,  being  influenced  by  torques,  is  called  precession. 
In  case  in  which  the  torques  come  from  the  friction  of  the  gimbal  bearings, 
then  the  angular  momentum  vector  again  precesses  ard.  thus  its  initial  orien- 
tation is  continuously  changing.  In  this  case,  we  speak  of  gyro  drift  which 
is  one  of  the  most  critical  errors  in  inertial  navigation. 


21 


Let  us  now  examine  how  three  gyros  can  either  stabilize  or  comand  a mo- 
ving platform  to  a desired  rate.  Consider  the  SDF-gyro  of  (Fig. (3b))  mounted 
on  a platform  which  generally  rotates  as  the  vehicle  moves.  Along  the  output 
axis  of  the  gyro  two  special  devices  are  set,  the  signal  generator (SG)  and 
the  torque  generator (TG) . The  former  has  the  ability  to  sense  torques  due  to 
the  rotation  of  the  gimbal  and  to  apply  them  as  restoration  torques.  The 
latter  can  generate  torques  according  to  the  rate  of  rotation  of  the  navi- 
gation frame.  For  example,  if  we  navigate  relative  to  an  earth  rotating 
frame,  then  the  torque  generator  is  mechanized  to  provide  torques  of,  appro- 
ximately, 15deg/h  and  the  signal  generator  senses  the  rotation  of  its  gimbal 
due  to  the  irregular  motion  of  the  vehicle  and  restores  it  in  order  to  iso- 
late the  platform  from  being  affected  by  vehicle  manoeuvres.  Ir.  case  we  na- 
vigate relative  to  an  inertially  non-rotating  reference  frame,  then  only 
the  signal  generator  operates  provided  that  the  gyro  frame  has  already  ocer. 
initially  aligned. 

Since  in  the  3D-space,  we  have  three  degrees  of  freedom,  it  is  understood 
tint  three  gimbals  provide  complete  isolation  of  the  platform’s  rotation 
witii  respect  to  the  turbulent  motion  of  the  moving  object.  Consequently, 
when  we  speak  of  three  on-board  gyves, we  mean  that  three  gimbals  preserve 
the  gyro  frame  to  change  its  attitude  as  the  vehicle  is  pitching,  rolling 
and  yawing. 

1.3.3  Jirror  model 

Various  gyre  error  models  have  been  suggested  in  an  effort  to  include  all 
possible  en-or  sources  in  a gyro  performance (e.g.  Britting,  1971,  p.74). 

The  use  of  a specific  model  depends  on  the  particular  problem  which  one 
faces.  As  we  have  already  explained,  our  gyros  command  the  platform  to 
the  earth’s  rate  and  thus,  the  error  model  should  be  selected  such  as  to 
includ’  the  instrumental  misfunctions  in  generating  the  earth’s  rotation. 

Let  us  suppose  that  we  have  on-board  the  moving  vehicle  three  SDF-gyrcs 
which  materilize  the  gyro  frame  generating  the  earth’  s rate.  Since  it  is 
technically  inpossible  to  construct  a pure  orthogonal  gyro  frame,  then  the 
resulting  frame  is  a non-orthogonal  or  quasi -orthogonal  frame.  Consequently, 
the  torques  are  applied  through  a non-orthogonal  frame  and  this  should  be 
taken  into  account.  Furthermore,  the  torques  are  experienced  through  the 
torque  generator  which  have  scale  factor  uncertainty.  Therefore,  the  signal 
for  the  earth’s  rotation  comes  erroneous  and  should  anyway  be  adjusted  be- 
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fore  it  is  applied,  taking  into  consideration  the  two  error  sources  just 
mentioned.  The  analysis  on  the  gyro  error  model  is  given  in  section  7.1, 


r 
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2.  The  fundamental  equation  of  inertial  navigation 

By  the  term  furriamental  equation  of  inertial  navigation,  we  mean  that 
equation  whose  solution  estimates  the  running  values  of  the  moving 
object.  In  our  analysis,  geocentric  coordinates  are  the  desired  output 
of  its  solution.  Consider  the  one  dimensional  navigation  exaiple  of  a 
train.  Its  fundamental  navigation  equation  is  distance  equals  velocity 
(considered  as  constant)  multiplied  by  time.  Measuring  time  by  an  on- 
board clock,  we  can  estimate  the  instantaneous  position  of  the  t'-ain 
with  respect  to  a convenient  initial  point. 

The  final  expression  of  the  fundamental  equation  depends  on  the  coordinate 
frame  in  which  it  is  coordinatized  as  well  as  the  motion  of  the  platform 
frame  relative  to  the  navigation  frame  used.  For  example,  if  the  navigation 
takes  place  ii:  an  inertial  coordinate  frame  and  the  platform  is  inertial ly 
stabilized,  then  the  navigation  equation  assumes  its  simplest  form.  In  ca- 
se in  which  the  platform  is  constrained  to  rotate,  then  the  navigation 
equation  includes  more  terms  such  as  accelerations  of  Coriolis  and  centri- 
petal type. 


We  shall  try  to  derive  the  fundamental  equation  of  inertial  navigation  ir. 
the  general  case  in  which  the  navigation  frame  rotates  inertially  and  the 
platform  frame  does  also  the  same.  Then,  it  is  easy  to  specialize  it  to 
any  desired  simplification.  We  finally  remark  that  the  navigation  equation 
is  valid  in  any  coordinate  system,  since  it  is  nothing  else  than  a vector 
equation.  It  comes  straightforward  from  the  total  derivative  of  a vector 
quantity.  Consider  two  coordinate  frames,  a fixed  one  and  a moving  E-,. 
The  transformation  between  them  can  be  represented 


£ 


2 


R £ 


1 


Differentiating,  we  get 


d£2  = dR  £j+  R dfi,  = dR  £1  = dR  RT  £2  = £2 

since  d£^  =0  being  fixed  and  n is  known  as  the  Cartan  matrix.  Now,  consider 
a vector  which  can  be  expressed  in  the  two  coordinate  frames  E^.and  E2  as 


a E1 
aiH 


bi  fi2 


where  a^  and  tn  the  coordinates  of  the  vector  A and  a^Ej  , b^  its  coirpo- 
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nents  in  E1  and  E2  frames  respectively.  Differentiating  again  we  get: 


(2.1) 


dX  ■ da^E^  + a^dE^  = db^  ^2  + dE2  or 
d 1 = da^  = dbt  E2  ♦ bt  Q E* 


Eq .(2.1)  can  be  interpreted  in  the  following  way.  The  derivative  of  a 
vector  with  respect  to  an  inertial  frame  is  equal  to  its  derivative  with 
respect  to  a moving  frame (such  as  one  fixed  at  the  earth's  centre)  plus 
the  rotation  of  the  moving  frame  with  respect  to  the  inertial  one  multi- 
plied by  the  vector  itself.  The  last  term  is  called  the  velocity  of  fol- 
lowing. We  state  again  this  important  conclusion  using  notation  valid  in 
what  follows  as: 


(2.2) 


dNr 


ar  * V 


dir 

it- 


where  I stands  for  the  inertial  frame  and  N for  the  navigation  one. Now 
let  us  derive  the  fundamental  equation  of  inertial  navigation  using  eq.(2.2). 
Consider  a moving  vehicle  0 connected  by  a vector  R to  the  navigation  fra- 
me N and  by  r to  the  inertial  frame  I (Fig. 4).  We  remind  again  that  R termi- 
nates at  the  platform's  mass  centre  and  when  we  speak  of  a moving  vehiclq, 
we  mean  only  a point,  its  platform’s  mass  centre.  From  the  vectors  definition, 


25 


where  p the  separation  of  the  origins  of  the  inertial  and  navigation 
frames.  Differentiating  with  respect  to  the  inertial  frame  we  get: 


(2.4) 


*i  ♦ “ W *4  * *n  " -ar  ?i 


(2.5)  \ * “ -3T"  K 


(2,6)  c SJ  )j  -ft  ♦ - c aj  )j  - * - aj  * ( ♦ aj A * ) • c *4 


Ec.  (2.7)  expresses  the  inertial  acceleration  of  the  platform’s  mass 
centre  as  a function  cf  the  acceleration  and  velocity  of  the  vehicle 
with  respect  to  the  inertial  one.  For  further  reference,  ej  r2.7)  is 
written  in  the  form 

(2.8)  = or  ♦ \ ~ ^ ^ “ R ♦ nj  M a*.  - t ) 


It  is  worth  noting  that  the  inertial  acceleration  is  the  sum  of  the 
gravitational  and  non-gravi rational  acceleration.  Consequently 

(2.9)  rj-X  + £ 

and  eq.(2.8j  assumes  now  the  form 


(2.10) 


1 S,  • a;  - % * 4J  * s£ 


If  we  want  to  estimate  geocentric  coordinates,  in  case  in  which  the  na- 
vigation frame  is  linked  to  the  earth,  then  we  have  to  integrate  eq.f2.10j 
doubly  with  respect  to  time.  But, as  one  could  foresee,  there  exists  some 
difficulties  in  performing  the  integration,  namely: 

a)  the  angular  velocity  of  the  navigation  frame  reffered  tc  the  inert id 
space  must  be  known.  Consequently,  it  is  intuitively  understood  that  t-e 
navigation  frame  should  be  linked  to  a body,  with  respect  to  which  we  like 
to  navigate,  which  has  known  angular  motion  characteristics  with  respect 
to  the  inertial  frame  of  reference.  Linking  the  navigation  frame  to  the 
earth’s  centre,  then  its  inertial ly  referenced  angular  velocity  is  a priori 
known  quantity  and  it  is  undoubtedly  a very  good  choice.  As  a matter  of 
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fact,  this  is  common  practice  in  inertial  navigation,  when  the  ravigation 
takes  place  near  the  earth’s  space. 

b)  the  inertially  referenced  angular  acceleration  of  the  navigation  frame 
must  be  known.  But,  in  view  of  the  obtained  accuracies  of  the  up-to-date 
navigation  systems,  this  quantity  is  very  small  and  could  be  neglected. 

c)  the  acceleration  of  the  navigation  frame  with  respect  to  the  inertial 
one  is  involved  in  eq.C2.10).  This  really  a problematic  quantity  as  far  as 
its  measurements  can  be  acconplished.  To  motivate  our  discussion,  consider 
the  centre  of  mass  of  the  universe.  This  point  satisfies  the  requirement 
for  constant  speed  rectilinear  motion,  that  is,  it  is  an  inertial  point. 
Now,  centering  three  axes  at  that  point  with  known  space  directions,  we 
get  an  inertial  frame  of  reference.  But  how  to  measure  the  acceleration 
of  the  navigation  frame  fixed,  say,  at  the  earth’s  centre  of  mass  with 
respect  to  the  inertial  frame?  And  even  if  we  continuously  approximate 
the  inertial  frame  with  other  quasi -inertial  ones  coming  closer  and  clo- 
ser to  the  navigation  frame,  certain  difficulties  in  measuring  the  iner- 
tial acceleration  exist  again.  To  overcome  this  cumbersome  situation,  let 
us  divide  the  universe  into  two  regions,  an  internal  region  and  an  exter- 
nal one.  Thus,  the  gravitational  acceleration  acting  on  the  moving  object 
is  the  sum  of  the  gravitational  acceleration  of  the  internal  and  external 
regions  or 


5 = 


If  we  consider  now  that  the  external  gravitational  field  Ggxt  is  uniform, 
that  is  to  say,  it  has  the  same  direction  and  magnitude  everywhere  in  the 

l»  t! 

external  region,  then  every  body  moving  in  the  external  region  will  get 
the  same  gravitational  acceleration  coming  from  the  external  region.  Since 
in  this  region  the  non-gravitational  forces  on  the  moving  body  are  zero, 
then  we  can  state 


In  view  of  this  result,  eq.(2.10)  reads 
*N  * * + Cint 

Now,  a very  reasonable  question  arises :with  what  criteria  one  could  divide 
the  universe  into  these  two  regions?  where  a*e  the  boundaries  of  them?  The 
answer  becomes  simple  in  case  we  know  in  what  kind  of  navigation  we  are 
involved.  In  our  studies,  the  navigation  takes  place  near  the  earth’s  sur- 
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face,  that  is  to  say,  we  are  interested  in  terrestrial  navigation.  Con- 
sequently, the  boundary  between  the  external  and  internal  regions  is  so- 
me kilometers  above  the  earth’s  surface  and  as  such  the  gravitational  field 
of  the  internal  region  represents  only  the  earth’s  gravity  field.  The  next 
to  the  mentioned  one  strong  gravity  fields,  namely  the  moon’s  and  sun’s 
fields,  are  neglected  being  quantities  up  to  the  order  of  10”7  the  earth’s 
gravity  field  intensity.  Finally,  the  fundamental  equation  of  inertial  na- 
vigation, in  case  of  terrestrial  navigation,  reads: 


(2.11) 


A + It, 


■♦I  "♦ 

24 ' \ 


Bi 

E 


where  X:  the  apparent  acceleration  of  the  moving  vehicle 
$p:  the  earth’s  gravity  field 

: the  Coriolis  acceleration 
nl  * (2  p 'Rj : the  centripetal  acceleration 

and  the  navigation  frame  is  centered  a:  the  earth’s  centre  of  mas:>. 


Eq .(2.11')  stands  for  what  we  are  after.  It  represents  the  fundamental  eqj- 
ation  of  inertial (terrestrial)  navigation  being  investigated  to  visual ice 
the  capabilities  of  a gradicmoter-aidcd  inertial  navigation  system  to  esti- 
mate instantaneous  geocentric  coordinates  of  a moving  vehicle  inside  the 
earth’s  space. 

It  is  practically  worthwhile  to  give  some  indicating  numbers  concerning  the 
magnitude  of  the  Coriolis  and  centripetal  acceleration.  Consider  f.  in 
eq. (2.11)  to  represent  the  earth’s  angular  velocity  and  an  aircraft  which 
moves  with  velocity  of  500  km/h.  Then,  the  Coriolis  acceleration  accounts 
for 


2G£  •*£  * 2 (7,29)0 O' S)rad/sec  50Ckm/h*  4-lO'3g 


where  g the  earth’s  gravity  field  intensity.  It  is  understood,  that  the  Co- 
riolis acceleration  magnitude  depends  upon  the  velocity  of  the  moving  t- 
hicle.  The  centripetal  acceleration  for  an  earth  bound  region  accounts  for 
almost  the  same  magnitude  given  above. 
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3,  Reference  frames  used  in  inertial  navigation 

The  discussion  about  any  inertial  navigation  configuration  intuitively 
includes  the  introduction  of  some  coordinate  frames  which  are  either 
ihstrumented  on-board  the  moving  vehicle  or  used  as  reference  frames 
and  as  such  the  system  output  is  related  to  them.  The  inclusion  and 
orientation  of  these  frames  depend  upon  the  particular  purpose  of  the 
navigation  system,  but  there  are  some  of  them  which  have  to  be  anyway 
implemented,  since  they  play  fundamental  role  in  the  whole  mathematical 
analysis. 

The  coordinate  frames  used  in  inertial  navigation  could  be  classified 
into  two  distinct  categories: 

3.1:the  external  coordinate  frames, which  are  linked  to  bodies  other  than 
the  moving  vehicle  and 

3. 2: the  internal  or  on-board  coordinate  frames,  which  are  linked  to  the 
object  being  navigated. 

In  the  first  category  two  particular  classes  can  be  defined: 

3.1.1: the  inertial  frame .where  the  term  stands  for  a number  of  quasi-iner- 
tial  frames  approximating  the  absolute  one.  Concerning  the  accuracy  of  the 
up-to-date  inertial  navigation  systems,  it  is  not  reasonable  to  desperately 
include  the  use  of  the  one  real  inertial  system.  Any  of  the  first  approxi- 
mations can  be  undoubtedly  considered  as  satisfactory.  One  could  adopt  a 
reference  system  with  origin  at  the  earth’s  centre  of  mass (included  its 
atmosphere)  and  with  the  following  axes  orientation: 

3-axis:  towards  the  instantaneous  rotation  axis  of  the  earth 

1- axis:  towards  the  vernal  equinox 

2- axis:  completes  the  right-handed  orthogonal  system 

Such  a quasi-inertial  frame  is  a common  choice  in  inertial  navigation  ap- 
plications, since  its  angular  velocity  with  respect  to  the  fixed  stars  ’S 
a negligible  small  quantity,  if  it  is  conpared  to  the  short  period  of  ope- 
ration of  a navigation  system.  The  frame  discussed  plays  the  role  of  the 
inertial  frame  in  the  present  work,  the  inertial  property  being  under- 
stood in  the  given  reasoning. 

3.1.2  : the  navigation  frame.lt  is  defined  as  that  frame  in  which  the  final 
output  of  the  navigation  system  is  coordinatized.  A lot  of  choices  could 
bring  a lot  of  navigation  frames  into  picture,  the  particular  choice  depen- 
ted  upon  the  purpose  of  the  navigation  system  under  consideration.  In  our 


analysis,  we  like  to  estimate  geocentric  coordinates  and  therefore  the 
navigation  frame  must  be  linked  to  the  earth.  The  orientation  of  that  frane 
is  as  follows: 

3-axis:  towards  the  rotation  axis  of  the  earth  considered  as  a solid 

body.  The  axis  is  fixed  at  the  time  instant  at  which  the  navigati- 
on starts.  If  we  like  to  refer  to  the  instantaneous  rotation  axis 
of  the  earth,  then  polar  motion  must  be  taken  into  account. 

1- axis:  towards  the  intersection  of  the  equator  fixed  when  the  naviga- 

tion begins  and  the  Greenwich  meridian  defined  by  a set  of 
world-wide  , well -distributed  astronomical  stations. 

2 - axis:  completes  the  right-handed  orthogonal  frame 

The  origin  of  the  navigation  frame  is  at  the  earth's  mass  centre  included 
its  atmosphere . 

Now,  let  us  discuss  the  internal  or  on-board  frames.  The  description  inclu- 
des all  or  at  least  most  of  the  or. -board  frames,  but  this  does  not  mean  that 
they  have  to  be  anyway  implemented. 

3.2.1  the  computation  frame.lt  exists  only  in  the  on-board  counter  memo- 
ry and  it  serves  as  that  frame  in  which  the  fin3l  manipulation  of  the  navi- 
gation system  equations  is  carried  out.  Having  decided  to  estimate  geocen- 
tric coordinates,  then  the  computation  frame  is  set  to  the  same  orientation 
as  the  navigation  frame  and  remains  parallel  to  it  everafter. 

3.2.2  the  mechanized  or  ideal  platform  frame. It  is  the  frame  with  respect 
to  which  the  gyro  and  accelerometer  frames  are  held  constant.  As  the  word 

It  II 

ideal  explains,  that  frame  is  practically  distorted  due  to  various  instru- 
mental reasons  e.g.  gyro  drift,  initial  misalignment  etc.  The  ideal  platform 
frame,  in  some  applications,  may  differ  from  the  navigation  frame,  the  moti- 
vation for  that  being  what  we  like  to  take  out  of  the  system.  In  our  analy- 
sis, the  ideal  platform  frame  is  identical  wi~h  the  computation  frame,  be- 
cause there  is  no  reason  to  assume  the  opposite. 

3.2.3  the  actual  platform  frame.  The  initial  platform  misalignment  cavses  the 
first  effect  and  as  the  time  passes  the  gyro  drift  drives  the  ideal  plat- 
form frame  away  from  its  desired  orientation.  The  resulting  frame  is  the 
actual  platform  one.  The  drifting  platform  affects,  finally,  the  orientation 
of  the  gyro  and  accelerometer  frames  which  lose  their  ability  to  ‘stabilize 
the  platform  and  sense  correctly  the  apparent  acceleration  components  respec- 
tively. 
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3.2.4  the  gyro  frame.lt  is  by  far  the  most  important  configuration  in 
any  navigation  system.  As  a matter  of  fact,  without  gyros  there  is  no 
way  to  navigate  a moving  vehicle (at  least  with  the  up-to-date  technolo- 
gy) . The  three  orthogonal  output  axes  of  the  on-board  gyros  construct 
the  sorcalled  gyro  frame  which  possesses  the  fundamental  principle  to 
preserve  its  orientation  regardless  the  turbulent  motion  of  the  carrier. 
Consequently,  every  distortion  of  the  platform’s  desired  orientation  can 
be  sensed  by  the  gyros  and  applied  as  restoration  torque  to  correct  the 
platform  to  its  proper  attitude.  Unfortunately,  the  gyro  frame  drifts  with 
time  and  this  complicates  the  navigation  analysis. 

3.2.5  the  accelerometer  frame. The  three  input  axes  of  the  on-board  acce- 
lerometers instrument  the  accelerometer  frame  along  which  the  apparent 
acceleration  is  sensed.  The  frame  is  affected  by  its  non-orthogonality 

as  well  as  the  drifting  platform. 

3.2.6  the  gradiometer  frame. As  in  the  gyros’  and  accelerometers’  case, 
the  axes  of  the  three  gradiometers  construct  the  gradiometer  frame.  It  is 
worth  noting,  that  this  frame  comes  into  picture  only  when  one  likes  to 
actually  measure  gravity  gradients.  Besides,  gradiometers  measure  pure 
gravity  gradients (without  inertial  disturbances)  if  and  only  if  they  fol- 
low an  inertial  orientation.  Consequently,  the  gradiometer  frame  must  be 
inertial ly  stabilized. 

Regarding  the  definitions  of  the  above  frames,  the  following  conments 
might  offer  clarity  and  understanding  in  what  follows: 
a)  The  aim  of  our  navigation  system  is  to  confute  the  instantaneous  geo- 
centric coordinates  of  the  moving  vehicle.  In  order  to  be  exact,  we  must 
declare  that  when  we  speak  of  a moving  vehicle  , it  must  be  understood  as 
only  the  centre  of  mass  of  the  moving  platform.  All  the  mathematical  rela- 
tionships as  well  as  the  estimated  geocentric  coordinates  refer  to  this 
point.  But,  there  are  two  arguments  against  this: 

I:  The  on-board  instmnentation  cannot,  of  course,  be  concentrated  on  the 
the  mass  centre  of  the  platform  and  it  is  actually  distributed  arotuid  this 
point  in  a radious  of  a few  meters. 

II:  The  gradiometers  are  tremendously  sensitive  instruments.  Thus,  it  is 
strongly  recommended  that  they  ought  to  operate  in  the  tail  of  an  aircraft, 
in  case  of  aircraft  navigation,  so  as  to  avoid  influence  of  the  movements 
of  the  personnel  etc. 

The  above  two  arguments  cause  no  trouble  because  the  centrifugal  accelera- 
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tions  created  from  the  distribution  of  the  instrumentation  around  the 
platform’s  centre  of  mass  (which  should  be  otherwise  induded)are  now 
dropped  of  the  navigation  system’s  equations  as  negligible  small  quan- 
tities. 

b)  It  is,  finally,  worth  repeating  that  the  orientation  of  the  navigation, 
confutation  and  ideal  platform  frame  is  the  same.  This  gives  simplicity, 
to  some  extent,  on  the  derivation  of  the  system’s  navigation  equation  and 
is  by  no  means  an  assumption. 
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irability  of 


radients  and  inertial  disturbances 


The  subject  of  this  analysis  is  of  principal  importance  in  case  in  which 
gradiometers  are  used  on-rboard  a moving  vehicle  to  measure  gravity  gra- 
dients. The  central  question  arisen  is  whether  we  can  measure  pure  gra- 
vity gradients  while  we  are  moving  or  not.  In  other  words,  we  have  to 
fiiri  out  what  is  the  output  of  a gradiometer  being  affected  by  the  tur- 
bulent motion  of  the  carrier. 

It  is  more  than  seventy  years  that  all  of  us  have  been  benefited  by  the 
Einstein’s  principle  of  equivalence  according  to  which  the  acceleration 
field  is  equivalent  to  the  gravitational  field.  It  is  worth  noting  here, 
that  this  principle  holds  only  approximately  and  locally  and  only  in 
that  sense  we  cannot  separate  gravitational  from  acceleration  effects 
(Fock,1964).  Let  us  therefore  examine  in  what  mode  a gradiometer  measu- 
res the  gravity  gradients  of  the  gravitational  field  in  which  it  operates. 

Consider  again  two  frames,  an  inertial  one  denoted  by  X and  an  accelerated 
with  respect  to  the  inertial  denoted  by  x.  The  relationship  between  then 
reads 


xi  • Vj  * Dt 

where  represents  the  relative  orientation  of  the  two  frames  and  the 
displacement  of  the  origin  of  the  moving  frame  with  respect  to  the  inertial 
origin.  Performing  in  exactly  the  same  mode  as  in  the  derivation  of  the 
fundamental  equation  of  inertial  navigation  and  representing  the  force  per 
unit  mass  by  F,  then  we  get 


dA.  ■ dx • 
-1&  IX1 


d2\j 
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Let  us  now  determine  the  terms  involved  in  eq.£4.1).  Using  the  Kronecker 
delta  definition  for  the  rotation  matrix  we  ootain 

Aik  Aij  ‘ 6kj 

and  by  differentiation 

dAik  dA.. 

(4.2)  —3^  • »ik  - 0 


^kj^jk 


33 


with  the  obvious  substitution  representing  the  instantaneous  inertiallr 
referenced  angular  velocity  of  the  moving  frame.  Differentiating 
with  respect  to  time  we  get 

,,  41  +A  Aj 

(4  -41  -nr — or  -at  *Ai  k 

and  according  to  eqs(4.2)  and  (4.3),  we  write 
dA  . dA  , 

<4-s>  hfyj  ■ l\‘  TT 

Consequently,  eq.(4.1)  is  now  written  as 


(4.6a) 


dx.  dn... 


where  f^  represents  the  coordinates  of  tne  gravity  vector.  E«-..;4.6a)  does 
not  assume  the  Newtonian  form  and  if  it  is  stretched  to  do  so,  wt.  obtain 


(4.6b) 


a a:- 

— t — ■ a,  » f.+x  r — x-f 

dtT  1 1 J V Ct  IK  JK 


where  the  ten  dx^/dt  has  been  dropped  out,  since  the  moving  instruments 
are  mounted  on  the  platform  and  therefore,  there  is  no  relative  motion. to. 
(4.od)  is  solved  witn  respect  tof.=  s\’/  and  then,  successive  space  dif- 
ferentiation yields  the  gravity  gradients  aid  the  third  gravitational  de- 
rivatives respectively.  Consequently, 

7'  3v-  ^ d2D, 

(4  ^ X'W-f*a:kVTr 


FV* j 


*Vv*T 


?A  o.'i  . 

Tc-  -(  St 


From  eqs (4 .7) , (4 .8)  and(4.S)  the  following  conclusions  can  be  drawn: 
a)  The  structural  difference  between  the  gravitational  and  acceleration 
field  allows  the  third  derivative  of  the  gravitational  field  to  be  free 
from  inertial  disturbances.  If  an  instilment  could  be  designed  to  measure 
the  third  derivatives,  then  its  measurements  would  be  pure  quantities  of 


f 
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the  gravitational  field,  the  motion  of  the  vehicle  having  no  effect  on 
them  whatsoever. 

b)  The  gravity  components  are  always  affected  by  inertial  disturbances. 
This  explains  the  fact  why  a gravimeter  cannot  perform  gravity  measure- 
ments while  it  is  moving. 

c)  When  the  platform  rotates  and  translates  in  space,  then  gravity  gra- 
dients are  mixed  with  inertial  disturbances.  The  way  out  is  to  inertially 
stabilize  the  gradiometer  frame  and  then  the  second  term  in  eq.(4.8)  re- 
presenting the  inertial  disturbances  is  cancelled.  But  as  we  have  designed 
our  platform,  it  trucks  an  earth  fixed  geocentric  frame  and  therefore,  it 
rotates  with  earth’s  rate.  To  reconcile  the  constraints  posed  on  the  gra- 
diometers,  a second  platform  is  set  on  the  moving  vehicle  which  is  iner- 
tially non-rotating  and  on  this  platform  the  gradiometers  perform  the  pu- 
re gravity  gradients  measurements. 

The  above  discussion  was  motivated  by  the  following  two  very  important 
reasons : 

1)  If  one  is  keen  on  asking  what  a gradiometer  measures,  then  eqs(4.7), 
(4.8)  and  (4.9)  give  the  straightforward  answer.  It  measures  what  one  li- 
kes as  well  as  what  one  dislikes.  It  is  in  one’s  choice, merely,  to  exclude 
inertial  effects  by  posing  the  constraints  which  come  so  clearly  from 
these  equations. 

2)  In  our  work,  we  like  to  estimate  geocentric  coordinates  referred  to 
an  earth-fixed,  non-inertial  frame.  But  our  discussion  has  proved  that 
gradiometers  must  measure  only  relative  to  an  inertial  frame  and  as  soon 
as  pure  gravity  gradients  have  been  obtained,  then  they  can  be  transformed 
to  any  other  reference  frame  such  as,  say,  to  the  earth-fixed  one.  Only 
under  these  lines,  it  is  well  understood  that  the  separability  of  gravi- 
tational and  inertial  effects  is  of  great  inportance,  when  gradiometers 
are  called  to  measure  only  pure  gravity  gradients  on-board  a moving  ve- 
hicle. 
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5.  The  system’s  driving  error  function 
5.1  General  considerations 

It  is  well-known  that  inertial  navigation  systems  are  burdened  with  a 
large  ammount  of  errors  which,  in  some  cases,  turns  on  to  be  intolerable. 
For  example,  the  vertical  channel  of  an  autonomous  navigation  system  is 
so  unstable  that  the  system  is  practically  out  of  its  reasonable  operation 
as  far  as  this  channel  is  concerned,  after  a short  time  of  operation.  The 
general  instability  of  an  autonomous  inertial  navigation  system  car.  be 
faced: 

a)  w~th  an  external  a^o:  various  instruments  have  been  suggested  and  used 
sc  far  in  an  effort  to  furnish  additional  on-board  measurements  with 
the  aim  of  reducing  a specific  channel  or  a whole  navigation  system’s 
instability.  For  tne  vertical  cliarnel  a barometer  or  an  altimeter 
could  accomplish  it  quite  succesfuliv  (Winter,  1974).  For  a whole  sy- 
stem, a velocitvmeter  or  a laser  equipment  which  measures  distances  or 
a gradiometer  or  a camera  taking  photographs  while  it  is  moving,  are 
some  of  the  external  aids  which  can  effectively  reduce  the  navigation 
system’s  errors  down  to  verv  reasonable  values. 

b)  with  a periodic  updating:  in  many  inertial  navigation  applications,  the 
system  becomes  rather  quickly  very  unstable  due  to  the  presence  of  big 
errors.  The  remedy  for  that  is  the  application  of  certain  mathematical 
methods, i.e.  Kalman  filtering,  to  update  the  system  at  a desired  observa- 
tion point  and  to  "start"  it  fresh  henceforth.  It  is  understood,  that 

t the  system’s  updating  takes  place  rather  often,  its  periodicity  being 
depended  upon  how  quick  the  errors  propagate,  the  sampling  interval  et^. 

Our  aniysis  deals  with  the  first  method  described  above.  Sc  far  the  earth’s 
gravity  field  is  approximated  by  the  gravity  field  of  a sphere  uc  an  ellip- 
soid of  revolution.  This  is,  of  course,  only  an  approx.ir.iat ion  of  the  reali- 
ty which  in  its  turn  contributes  errors  to  u.~-  nrvigatior.  system  peril 
In  each  science  special  assumption s are  always  set  in  an  effort  tc  approach 
the  reality  as  closer  as  possible.  But,  in  view  of  the  hopeful  results  co- 
ming from  the  gradiometer  feasibility  studies  and  tests,  we  could  take  full 
advantage  of  them.  Now,  we  have  an  instrument  which  has  the  ability  to  me*- 
sure  the  gravity  field.  Thus,  why  to  make  such  a vital  assumption?  Finally, 
we  enrich  our  platform  instrumentation  by  setting  a number  of  gradiometers 
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to  measure  the  earth’s  gravity  field.  Certain  assumptions  have  of  cource 
to  be  made  concerning  gradiometer  accuracy,  values  of  the  gravity  gradi- 
ents in  space,  location  of  the  platform’s  instrumentation  etc. 

*5.2  The  simulated  navigation  equation 

Our  study  examines  one  of  the  simple  cases  of  an  inertial  navigation  sy- 
stem, trying  to  firxl  the  function  which  drives  the  error  budget  of  the 
system  under  consideration.  In  order  to  bring  the  fundamental  equation 
of  inertial  navigation  into  an  easily  simulated  form,  certain  assumptions 
have  to  be  made,  e.g.  absence  of  instrumental  errors.  We  believe  that  ma- 
king these  assumptions,  we  do  not  overshadow  problems,  but  we  present  a 
simplified  example  to  draw  a very  interesting  conclusion,  namely,  which 
is  the  worst  contributed  factor  in  an  inertial  navigation  system. 

Let  us,  now,  examine  in  detail  our  navigation  system.  The  moving  inertial 
platform  is  on-board  an  aircraft  and  is  composed  of  the  following  instru- 
mentation: 

a)  three  orthogonally  mounted  single-axis  accelerometers  in  order  to  mea- 
sure the  components  of  the  apparent  acceleration  vector  resolved  into 
their  axes. 

b)  three  orthogonally  mounted  gyroscopes  which  can  isolate  the  platform 
from  the  turbulent  motion  of  the  aircraft  and  preserve  the  orientation 
of  the  platform  frames  with  respect  to  the  inertial  frame. 

c)  three  orthogonally  mounted  spherical  gradiometers  which  furnish  the  mea- 
surements of  the  entire  gravity  gradient  tensor. 

Taking  into  account  that  the  moving  platform  does  not  rotate  with  respect 
to  the  inertial  frame,  which  is  in  this  case  the  navigation  frame,  then 
all  terms  in  the  general  navigation  equation (eq. (2. 1 1))  containing  the 
angular  velocity  term,  are  simply  dropped  out.  Consequently,  the  navigation 
equation  covering  our  system  assumes  the  simple  form 

(5.2.1)  R1-A*GI 

where  R denotes  the  instantaneous  inertial  position  of  the  aircraft,  sub- 
script I shows  the  inertial  reference,  A is  the  apparent  acceleration  and 
G the  gravity  conponents.  For  simplicity,  the  inertial  subscript  is  dropped 
henceforth.  Eq. (5.2.1)  can  be  written  in  component  form  as 
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x‘VSc 

(5.2.2)  Y-yGy 

Z«A  ♦G 
Z Z 

where  X,Y  and  Z the  inertial  coordinates  of  the  aircraft.  Since  we  do  not 
measure  gravity  components  but  gravity  gradients  and  the  inertiai  accele- 
ration terms  call  for  undesired  integration,  eqs(S.2.2)  are  approximated 
in  the  following  reasoning: 

a)  the  inertiai  acceleration  components  could  be  substituted  by  the  inertial 
coordinates,  using  the  well-known  Stirling’s  approximation  formula  appli- 
ed  for  three  successive  points (Scheid,  196SJ 


(5.2.3) 


V 


yi-w, 

771 


b)  the  gravity  components  are  to  be  approximated  by  the  gravity  gradients 
using  a coordinate-free  derivation  which  is  valid (as  the  name  explains) 
in  any  coordinate  frame: 

(5.2.4)  Gi  « Gl_.*gradC1_;(Ri-Ri_1)*  higher  order  terms 
where  the  higher  order  terns  are  neglected. 

Consequently,  eq.(S.2.2)  in  view  of  the  approximations  made  could  be  re- 
written in  the  form: 


<xi ♦r2VW*t‘2  ■^/%+CXXi./Xi-Xi-t)+GAC.1^i-Yi-l)+GXZ1./2i-2i 

quantization  error  ‘ 

(5.2.5)  Of^-ZY^ytt-2  ■\.+Si+Grxi./xi-xi-i)+GY3'i./Yi'Yi-i]+Gi7i..r::i-:i-i)+ 

quantization  error 

quantization  error' 

Qiantization  error  studies  are  performed  in  section  6.  Eqs(5.2.6)  show 
a rather  peculiar  phenomenon.  When  the  moving  platform  is  at  the  point  i 

*i*i  ■ ^2^^i./vxiM)%i/iyi)tGxii./^-Vi))+2xi-xi-i4  q-e- 

($.2.6)  Yi4,  - At2(Ayi^i^i_i(Xi-Xi.1)^i_i(Yi-Yi.1)^Zi_i(Zi-Z..1))*2Y.-YI.r  q.e. 


(Z.-Z.  *)♦ 

14,4 
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Vi  " At2(Azi*GZi+GZXi_1(Xi_Xi-1J^GZYi.1(Yi‘Yi-1'1+GZZi_1(Zi'Zi-l^+2Zi‘Zi-i+  q,e- 

and  the  on-board  instrumentation  measures  apparent  acceleration  and  gra- 
vity gradients (as  a matter  of  fact  the  measurements  are  performed  oetween 
the  points  (i-1)  and  i) , then  the  system  estimates  the  inertial  coordina- 
tes of  the  (i+1) -point.  Therefore,  the  platform  behaves  like  a "moving  win- 
dow” and  there  are  two  general  ways  one  could  face  this  situation: 

a)  the  Boundary  Value  Problem  of  Inertial  Navigation: given  the  inertial 
coordinates  of  the  first  and  the  final  points,  then  the  system  can  esti- 
mate all  other  points,  but  only  off-flight,  that  is,  post  mission.  In  so- 
me applications,  this  has  certain  advantages.  For  example,  consider  a pho- 
togrammetric  airplane  taking  photographs  of  an  ar^z.  During  the  flight,  the 
instantaneous  coordinates  of  the  moving  vehicle  are  not  of  direct  interest. 
But,  after  the  mission  has  been  accomplished,  the  desired  photo-manipulation 
needs  the  coordinates  of  the  camera  from  where  all  photographs  have  been 
taken.  Thus,  application  of  eqs (5.2.6)  simultaneously  for  all  in  flight 
points  gives  the  required  values. 

b)  the  Initial  Value  Problem  of  Inertial  Navi  gat  ion -.when  the  first  two 
points  of  the  navigation  path  have  known  coordinates,  then  eqs (5. 2. 6) 
estimate  the  inertial  coordinates  of  all  other  points  to  come  in  flight. 

In  many  practical  inertial  navigation  applications,  as  in  terrestrial  na- 
vigation, we  are  most  interested  in  knowing  at  any  time  instant  where  we 
are  and  therefore  this  case  suits  to  our  analysis.  Schematically,  our  na- 
vigation problem  is  illustrated  in  Fig. (5). 
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Fig. (5):  The  initial  value  problem  inertial  navigation 


In  Fig.  (5),  quantities  above  the  navigation  path  ar^to^be  measured  and 
under  it  are  known.  The  quantities  inside  the  brackets  underthe~naviga- 
tion  path  indicate  what  is  estimated  when  the  platform  is  at  that  point. 

Before  we  compute  the  matrices  involved  in  the  variance -covariance  expres- 
sions for  the  estimated  inertial  coordinates,  it  is  worthwhile  to  examine 
the  inpact  of  the  Laplace  condition  on  the  variances  of  the  measured  gra- 
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vity  gradients.  It  is  well-known  that  the  in-line  gravity  gradients  are 
connected  with  the  following  condition: 

(5-2.7)  Gz,  - •gxx‘(Y l 

known  as  the  Laplace  condition.  For  instructive  purposes,  we  rewrite  the 
gravity  gradient  tensor 

Sex  Scy 

Srx  j Sr  S z 

GZX  CZY  GZZ 

and  note  the  solenoidal  and  irrotaticnal  structure  of  the  gravity  field 
(Gf^Gji).  What  is  above  the  dotted  line  is  observed  and  certain  accura- 
cies are  later  to  be  assumed  for  them.  But,  the  vertical  in-line  gravity 
gradient (G^v)  which  enters  eq.  (5.2.7)  has  to  get  a variance  dependent  on 
the  other  two  gravity  gradients  variances.  Consequently,  trying  to  find 
the  relationship  between  the  in-line  gravity  gradients  variances,  we  write 

Gxx“Sa*°Sr 


SY"0G)«+Sr 

Gzz'"Gxx'Sr 

or  in  matrix  form 
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If  eq.  (5.2.8)  is  abstractly  written  as  Y=AX  and  the  law  of  propagation  of 
errors  is  applied,  then  we  obtain 
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(5.2.9) 
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"Vyy^ 


-<3gXx"°CVyGXX  ~<'GXXGYY  7gYY  0gXX  0(Vy  aGXXGrY| 


Equating  the  respective  matrix  elements,  we  get  besides  the  obvious  iden- 
tities the  following  interesting  expressions: 

GzzGxx  °xx  SxSt 


(5.2.10)  °^2ZG*YY  0cVy  °GXXGYY 


oi  * ai  _ 

Gzz  Wy  ‘’xxWy 

The  above  equations  are  of  value  for  our  analysis  because 

a) when  we  consider  the  observables  to  be  uncorrelated,  then  we  cannot  as- 
sume the  same  for  the  vertical  gravity  gradient  and  thus  its  correlation 
with  respect  to  the  other  two  in-line  gradients  ought  to  be  considered. 

b) when  we  simulate  the  navigation  system,  certain  variances  have  to  be 
chosen  for  the  G^,  C^.  gradients.  Gz?  will  then  get  a variance  as  eq. 
(5.2.10c)  indicates. 

What  we  intend  to  examine  now  is  the  contribution  of  each  observable  of 
the  navigation  system  to  its  error  budget.  The  simulation  takes  place 
when  the  platform  is  at  the  2 -point  and  so,  the  equations  can  analytically 
be  presented  as: 

X3* t2  (AX,  2+GX1-‘GXX,  2{X2'X1  > *CXY, 2 CY2’Y1^GXZ1 z(Z2'Z1  »+2VX1 
(5.2.11)  Y3*t2(Ay  ♦Gy  *Sxi,CX2‘X1)+GYY12(Y2'Y1)+GYZ12(VZ1))’2VY1 

12  1 1 fc  * 


Z3.^(aZi^Gz^gZxi2(Vxi)-gzYi2(Y2-Yi)^Z2i2(Z2-Zi))42Z2-Z1 


In  the  system  of  eqs(5.2.11)  we  have  19  parameters,  namely 
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a)  At: the  time  interval  between  two  successive  observations (1) 

b)  : the  components  of  the  apparent  acceleration (3) 

c)  : the  gravity  components (3) 

d)  G^:three  in-line  and  three  off-line  gravity  gradients(6) 

e)  (X,Y,Z) 1 , (X,Y,Z)2:the  inertial  coordinates  of  the  first  two  points (6) 
Generally,  the  variance-covariance  matrix  of  the  estimated  inertial  coor- 
dinates can  be  written 

(S.2.12)  t y - G1^ 

where  r y indicates  the  variance-covariance  matrix  of  the  coordinates, 

Z^  the  variance-covariance  matrix  of  the  19  parameters,  G the  matrix  of 
the  derivatives  of  eqs(5.2. 1 1)with  respect  to  the  19  parameters  and  T 
denotes  the  transpose  matrix.  In  view  of  eqs (5.2.11),  the  simulated  sy- 
stem is  given  in  table  (1).  The  derivatives  of  eqs (5.2. 11)  with  respect 
to  the  parameters  involved  read: 

aZT"  “iit(AX12*S1*GXX12(X2'X1)+ScY12(Y2‘Y1)+GX212{VZ1^ 
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(5.2.14) 
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5.3  Results 

Eq. (5.2. 13)  is  our  simulated  navigation  equation  combined  with  eqs(5.2.14)‘. 
In  order  to  obtain  the  desired  results,  certain  accuracies  have  been  assu- 
med for  the  19  parameters  of  the  system.  The  simulation  covers  the  case  in 
which  the  observables  are  not  correlated.  When  an  observable  assumes  a ran- 
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ge  of  values,  the  rest  18  are  held  constant  so  as  the  influence  of  the 
observable  under  consideration  on  the  variances  of  the  inertial  coordi- 
nates could  be  deduced.  The  used  ground  values  for  the  19  parameters  of 
the  system  are  listed  in  Table  (2). 


Parameter 

Value 

Variance 

n 

0.1  sec 

0.001  sec2 

II 

0.1m/secZ 

0.01m2/sec4 

\ 

0.1 

0.01 

0.1 

0.5 

0.01 

0.5 

0.01 

% 

9.0 

0.01 

Sa. 

f 

-1500.E 

1.E 

Gxy) 

-0.2 

1. 

-15. 

1. 

s 

-1S00. 

7. 

S 

♦0.05 

1. 

G2Z 

♦3000. 

1. 

X1 

6200000  m 

I.m2 

Y, 

6000000 

1. 

Z1 

6500000 

1. 

x2 

6200030 

’•  i 

Y2 

6000030 

1. 

Z2 

6500030 

1. 

Table  (2) : Ground  values  for  the  parameters  used  in  the  simulation. 
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As  far  as  these  values  are  concerned  the  following  remarks  are  indi- 
spensable: 

a)  the  operation  of  the  platform’s  instrumentation  every  AtO.lsec. 
could  be  generally  considered  as  rapid,  but  actual  tests  on  aided  in- 
ertial navigation  systems  recommend  that  value (Denhard,  1977). 

b) the  values  for  the  initial  inertial  coordinates  are  selecteu  to  be  out- 
side the  earth’s  surface  having  30m  distance  apart.  An  aircraft  velo- 
city of  lOOOkm/h  is  assumed  for  the  moving  platform.  But,  generally 
speaking,  the  given  coordinate  values  have  not  any  special  meaning 
whatsoever. 

c) the  ground  values  for  the  gravity  gradients  have  been  taker,  from  (Tra- 
geser,  1971)  and  compared  to  those  given  by  (Reed,  1973).  It  is  worth 
noting  here  that  for  our  simulation  studies  only  indicated  values  for 
the  gravity  gradients  fulfill  the  requirements  of  the  analysis,  since 
we  do  not  derive  coordinates,  but  the  influence  of  each  observable  on 
the  variances  of  the  coordinates. 

For  the  computational  analysis,  a special  programme  was  written  in  the 
sense  of  eqs(5.2.12)  and  compiled  in  the  Siemens  Computer  at  the  Univer- 
sity FAF  Munich.  The  results  of  the  carried  out  simulation  arc  given 
below . 


8 


var . -cov . : 


4.0340  0.0040 

5.0240 


4.0740  0.0040 

S.0240 


4.S240 


6.0240 


9.0240 


0.0040 

5.0240 


0.0040 

5.0240 


0.0040 

5.0240 


0.0400 

0.04J0 

5.4201 


0.0400 

0.0400 

5.4201 


vir . -cov . : 

(m2) 

1 .0640 

0.0040 

0.0400 

5.0240 

0.0400 

5.4201 

1.2240  0.0040 

5.0240 


1.4240  0.0040 

S.0240 


C.0400 

0.0400 

5.4201 


0.0400 

0.0400 

5.4201 


0.0400 

0.0400 

5.4201 


3.0240 


0.0040 

5.0240 


0.0400 

0.0400 

5.4201 


0.0400 

0.0400 

5.4201 


0.0400 

0.0400 

S.4201 


0.0040 

0.0400  i 

5.0240 

0.0400  1 

21 .02*0  0.0040  0.0400  | 

5.0240  0.0400  ! 


var. -cov. : 

On2) 

5.0241 

0.0040 

0.0400 

5.0241 

0.0400 

5.4202 

5.0240 

0.0040 

0.0400 

5.0245 

0.0402 

S.4207 

At  (sec 


5.0240 

0.0010 

5.0240 

0.0400 

0.0400 

5.4201 

5.2200 

0.2000 

1.9993 

5.2200 

1.9998 

25.0213 

5.4199 

0.3999 

3.9995 

5.4200 

3.9998 

45.0225 

7.0194 

1.9996 

19.9976 

7.0197 

19.9991 

205.0321 

0.0400 

0.0450 

S.4326 


var . -cov . : 


5.0240  0.0040 

5.0240 


5.0240  0.0040 

5.0240 


5.0240  0.0010 

5.0240 


0.0400 

0.0400 

4.4301 


0.0400 

0.0400 

4.4701 


0.0400  l 
0.040G  i 


var. -cov. : {m  ) 


5.0240  0.0040 

5.0240  0.0400 


5.0240 


5.0240 


0.0040 
5. 


0.0010 

5.0240 


0.0100 


1.6201  I 


0.04X 

0.0400 


.8 


! S.00 


5.0240  0.0040  0.0400 

S. 


5.0240 


5.0240 


5.0240 


0.0040 

5.0240 


0.0040 

5.0240 


0.0400 

0.0400 

6.4201 


0.0400 

0.04X 

9.4201 


0.0040 

5.0240 


0.0400  , 
0.0400  I 
3.420'  ! 


j 

! 5.0240  0.0040  0 

5.024  0 0 

9 

.0400 

.0400 

.4201 

1 

1 5.0240 

1 

0.0010 

0.0400 

5.0240 

0.0400 

X (m  /sec 


var. 

, -cov . : 

oo 

5.0211 

0.0020 

0.0201 

5.0240 

0.0400 

5.4201 

5.0211 

0.0021 

O.C206 

5.0240 

0.0400 

5.4201 

5.0212 

0.0022 

0.0220 

5.0240 

0.0400 

5.4201 

5.021S 

0.0024 

0.0240 

5.0240 

0.0400 

5.4201 

i’“i-1.2  fm 


0.0740 

0.0040 

0.0*40 

0.5240 

0.0040 

0.0400 

| 

0.5240 

0.0400' 

0.0040 


.9 


0.0400! 


10.0240  0.0040  244.0042 

10.0240  244.0225 

4882.5818 
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From  the  obtained  values  the  following  conclusions  can  be  drawn: 

a) as  the  accuracy  of  the  clock-timer,  which  signals  the  time  interval 
At  and  alerts  the  platform  instrumentation  to  operate,  decreases,  then 
the  variance  of  the  Z-coordinate  gets  very  big  values.  Generally  spea- 
king, it  is  reasonable  to  assume  that  the  clock-timer  possesses  no  error 
in  determining  At(observe  also  the  increasing  values  of  the  covariances!). 

b) the  accuracy  with  which  the  inertial  coordinates  of  the  second  point 
are  known  causes  worse  instability  than  that  of  the  first  point.  The  va- 
riance of  the  X-coordinate  has  no  effect  on  the  derived  Z-coordinate  va- 
riance and  the  deviations  in  Z cause  bigger  instability  to  Z-variance 
than  that  due  to  X or  Y variations. 

c) the  accuracy  of  the  gravity  gradients  has  almost  no  effect  on  the  esti- 
mation of  the  inertial  coordinates.  This  is  really  something  surprising 
and  not  known  so  far.  It  can  be  explained  as  follows: if  one  performs  the 
tedious  manipulation  on  the  simulated  navigation  eq. (5.2.13),  one  could 
see  that  all  terms  containing  gravity  gradients  are  multiplied  by  the 
fourth  power  of  the  time  interval  At  and  by  other  small  quantities.  Thus, 
their  influence  is  strongly  reduced. 

d) the  variance  of  the  gravity  components  of  the  first  point  has  almost  no 
influence  on  the  variance  of  the  coordinates  of  the  third  point.  Even  if 
the  variance  is  equal  to  the  value  of  the  gravity  component ( ! ) , then  the 
variance  of  the  third  point  is  burdened  by  half  a meter  more. 

e)  acceleration  seems  to  have  no  great  influence  on  the  derived  accuracy 
of  the  third  point.  But,  as  we  shall  see  later,  acceleration  variation 
creates  the  biggest  errors  in  the  system  among  the  10  observables. 

So  far  the  analysis  has  given  the  first  results,  namely,  the  behaviour 
of  the  navigation  system  with  respect  to  the  parameters.  Next  we  try  to 
investigate  how  the  system  performs  going  from  point  to  point.  During 
that  process  the  observables  are  undergone  small  changes  except  for  the 
variances  of  the  coordinates  of  the  first  two  points  which  are  assumed 
to  be  zero.  The  same  computer  procedure  is  used  and  applied  for  some 
first  points  due  to  reasons  which  will  be  explained  later  on.  The  space 
traverse  gives  the  following  results: (see  next  page) 

5.4  Discussion 


First,  it  has  to  be  noted  that  the  statistical  analysis  up  to  now  has  not 
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Point 

mm 

era 

2 

°2 

3 

0.000016 

0.000016 

0.003314 

4 

0.000080 

0.000080 

0.016570 

S 

0.000352 

0.000352 

0.072908 

6 

0.001504 

0.001 504 

0.3115 

7 

O.OOt384 

0.006384 

1.3222 

8 

0.02706 

0.02706 

5.6036 

9 

0.1146 

0.1146 

23.7399 

10 

0.4855 

0.4855 

100.5665 

2.0566 

2.0566 

426.0092 

12 

8.7095 

8.7095 

1804.6066 



36.8966 

36.8966 

j 

7644.4123 

Table  (3):  A space  traverse  of  a gradiometer-aided  inertial  navigation 
system. 

included  any  error  model  of  the  platform  instrumentation.  But,  certain 
errors  do  exist  which  sometimes  turn  to  be  important.  Gyro  drift,  ini- 
tial platform  misalignment,  gyro  ner.-ortnogonality  etc.  are  some  errors 
which  in  a rigorous  statistical  analysis  have  to  be  modelled  and  taken 
into  consideration.  All  these  problems  are  thoroughly  examined  in  simu- 
lation II . 

The  performed  simulation  studies  did  have  the  objective  to  be  simple  and 
as  such  the  following  conclusions  can  be  drawn: 

a) first  of  all  it  is  clear  that  the  gradiometer-aided  inertial  navigation 
system  we  analyse  does  not  perform  well  as  far  as  the  accuracy  of  the  de- 
rived coordinates  is  concerned.  The  first  few  points,  say  up  to  the  tenth, 
could  be  reached  with  satisfactory  precision^see  Table(31).  Then  the  accu- 
mulation of  the  system’s  errors  becomes  so  high  that  the  performs.' 'e  of  it 
can  be  considered  as  unreasonable. 

bjraanipulating  the  expressions  of  the  variances-covariances  of  the  simu- 
lated equations (Table(1 )) , we  find  that  besides  all  other  terms,  there 
are  two  of  great  interest 

1)the  variance  of  the  (i-1) -point  multiplied  by  four  and 


2) the  variance  of  the  (i -2) -point 

in  case  the  inertial  instrumentation  is  at  the  i-point.  It  is  worth  noting 
that  both  terms  have  positive  sign. 
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c) the  answer  to  the  bad  behaviour  of  the  system  lies  on  the  conclusion 
b) . If  we  assume  that  the  first  two  sets  of  the  inertial  coordinates 
have  zero  variances,  then  all  other  terms  of  the  variance-covariance 
expressions  are  summed  to  give  the  variances  of  the  cbordinates  of 
the  third  point.  Four  times  this  variance  plus  the  variance  of  the  se- 
cond point  (plus  other  small  terms)  will  give  us  the  variance  for  the 
fourth  point.  The  same  procedure  is  applied  through  all  other  points. 
It  is  therefore  clearly  understood  that  the  following  approximate  law 

i-point  variance  = 4* (i-l)-point  variance+(i-2) -point  variance 

is  valid  for  the  space  traverse  under  consideration. Consequently,  an 
appreciable  percentage  of  the  navigation  system’s  error  budget  comes 
from  the  approximation  used. 

d) next  we  try  to  see  the  error  contribution  of  the  acceleration  and 
gravity  gradients  measurements  into  the  system.  To  that  objective  two 
additional  space  traverses  are  performed.  In  the  first  one  all  other 
quantities  except  the  apparent  acceleration  components  assume  zero 
variances  throughout  the  navigation  path,  Thus,  the  error  contribution 
of  the  function  under  consideration  can  be  deduced  and  compared  to 
the  system’s  error  budget  at  a selected  point.  The  same  procedure  is 
followed  for  the  gravity  gradients  case.  The  results  are  listed  below: 


Accelerometer 

Gradiometer 

Point 

Q 

X rvj 
it 

Q 

~<  ru 

w 

2 , 2, 
o2:(m  ) 

T2 

°X*°Y 

3 

o.to 

10'5 

2S0 

i0-10 

4 

0.50 

10'5 

1 .35 

10"9 

5 

0.22 

10*4 

5.94 

O 

< 

VO 

6 

0.94 

to*4 

2.53 

to-8 

7 

0.40 

to’3 

1.07 

to-7 

8 

0.T7 

N 

1 

o 

4.56 

to'7 

9 

0.72 

10*2 

1.93 

to"6 

to 

0.31 

to"1 

8.19 



to'6 

From  the  listed  results  it  is  evident  that 

1) the  accelerometer  error  contribution  into  the  system’s  error  budget  is 
about  7%  of  the  total  and 

2)  the  gradiometers  work  perfectly  well  in  such  a navigation  system.  Their 
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error  contribution  is  really  minimal. 

In  the  next  pages  we  give  some  representative  nomographs  which  illu- 
strate the  behaviour  of  the  parameters  of  the  system  relevant  to  the 
obtained  accuracy  of  the  inertial  coordinates. 

In  section  7 we  will  come  again  to  discuss  the  same  topic  but  this  time 
detailed  error  models  for  the  accelerometer  and  gravity  gradients  functions 
will  be  included  to  see  if  the  system’s  results  could  be  improved. 

Exptanatcay  note:  The  noltoieing  iym.boli*m  -Li  u*cd  in  the.  attacked  no- 

mogaaphi  e.g.  cv  (A..  ).  1 1 mean. i that  of  unde’igoa 

the  change i plctuaed  by  the  ae-ipective  Lcne  only  token 

Ay  a&4umc>6  a \anqe  o{,  value*. 

A12 


2.0 


O'O 


6. 


Inertial  measurement  unit  quantization  error  analysis 


6.1  The  general  equation 

Quantization  errors  constitute  a special  error  category  in  every  instru- 
mental package  and  must  be  clearly  distinguished  from  round-off  errors, 
commutation  errors  etc.  Quantization  errors  arise  from  the  fact  that  the 
inertial  measurement  units  measure  continuous  physical  quantities  such 
as,  say,  the  apparent  acceleration  vector  components.  Assuming  that  at 
each  observation  point  at  which  the  inertial  instruments  are  read  a quan- 
tization error  exists,  then  its  additive  net  effect  will  burden  the  de- 
sired output  of  the  navigation  system  with  an  additional  error.  In  our 
analysis,  the  system  computes  geocentric  coordinates  making  use  of  acce- 
lerometer and  gradiometer  measurements.  Expressing  position  with  respect 
to  the  measurables  and  attaching  to  them  different  quantization  errors, 
then  the  position  error  could  be  computed  treating  the  quantization 
errors  as  stochastic  quantities.  Regarding  the  up-to-date  accuracy  of 
the  inertial  navigation  systems,  it  will  be  seen  that  the  quantization- 
induced  position  error  is  not  oi  great  inportance (due  to  its  sm.il  1 mag- 
nitude). But,  there  is  one  case  in  which  that  error  could  be  of  great  inte- 
rest. As  we  have  already  seen,  all  inertial  navigation  systems  without 
exception  must  be  initially  aligned  to  a desirable  coordinate  system 
prior  to  their  mission.  In  some  cases,  the  initial  alignment  procedure 
can  consume  an  appreciable  long  time  (especially  in  commercial  flights) 
and  during  that  time  the  inertial  instrumentation  is  continuously  ope- 
rating and,  of  course,  gathering  quantization  errors.  These  errors  could 
cause  an  initial  misalignment  error  which  is  carried  through  the  entire 
mission. 

Consequently,  it  is  cur  belief  that  quantization  errors  must  be  always 
considered  to  preestimate  how  badly  they  burden  the  system’s  output.  In 
what  follows,  we  present  the  first  unique  grzliometer-aided  inertiai  na- 
vigation systems  quantization  error  analysis  using  the  fundamental  equa- 
tion of  inertial  navigation  used  in  the  preceding  simulation. 

As  we  have  already  seen,  the  formula  which  estimates  the  infertial  coordi- 
nates is  written  in  the  form: 
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where  A and  G the  non-gravitational  and  gravitational  acceleration  re- 
spectively, P indicates  position (in  eq.  (5.2.1)  R is  used  instead  of  P) 
and  n represents  the  nth-observation  point.  Stirling’s  formula  approxi- 
mates the  inertial  acceleration  as: 

. (6<1-2)  Pn  -at'2<Vr2P„+  pn-1> 

Eq. (6.1 .2)  solved  for  the  nth-inertial  coordinate  gives 

(6.1.3)  V2Pn-rPn-2+6tVit2Gn 

Our  policy  is  row  to  express  the  above  equation  as  a function  of  the 
first  two  inertial  coordinates  which  are  then  to  be  considered  as  error- 
less quantities.  In  order  to  find  the  general  formula  which  gives  the 
coordinates  of  the  n^-point,  we  write 

P.-2P,-P,-At2A,*At2G, 

(6.1.4)  J i 1 " l 

P4-3P2-2P)*At2(2A2+Aj)+At2(2G2*GJ) 

Consequently,  we  get 

Pn«  (n-1  )P,-  (n-2)P1+at2((n-2)A2+  (n-3)Aj+  • • (u-(n- 1 ) )+At2(  (n-2)G,+  in-5)G-+  • ■ 

(6.1.S)  ♦(n-(n-1))Gn_1) 

and  in  suimation  form 

n-l  . 

(6.1.6)  Pn-(n-1)P2-(n-2)P1*E 

m*2 

Now,  employing  the  free -coordinate  approximation  for  the  gravity  gradi- 
ents we  write 

(6.1.7)  G -G  . *gTadG  , (P  -P  ,) 

«i  m~ 1 m- 1 v m m- 1 1 

Expressing  the  above  equation  in  terms  of  the  gravity  components  of  the 
first  point  plus  other  gravity  gradients,  we  get 


i 

i 


l 
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•(6.1.8)  G^CgradC^^-P^) 


Consequently,  eq.  (6.1.6)  is  now  written 


(6.1.9)  Pn-(n-1)P2-(n-2)P1^(n-m){Ah|+  G,*"  gradG^,  CPk-Pk_1))  At2 


Assuming  now  that  each  time  the  accelerometers  and  gradiometers  are  read 
a quantization  error  is  committed,  we  express 


n-1 


(6.1 .10)  Fn^Pn.(n-l)P2-(n-2)P-^(n-m)  (Am+  )At2+G,?r  (n-nD^fgradG^,  CPjc-F-._1  )*q«_- )tt2 


m-2  k= 


where  rquantization-induced  position  error 


qm  '.accelerometer  quantization  error 


q^_1 ^radiometer  quantization  error 


Taking  into  consideration  only  the  impact  of  the  quantization  errors, 
then  the  position  error  assumes  the  form 


7 n_  1 


, n- 1 m 


(6.1.11)  i>P  ‘A:*-  l (n-n)o  ♦ At*  f (n-irj  £ q? 

n m-2  ™ m=2  • k«2  ‘ 


We  shall  derive  now  the  general  relation  for  the  variance -covariance  of 


the  position  error  6Pn  . Generally,  we  can  write 


(6.1.12)  cov,SPn  -F(6PJheP  -E(6P  )} 

n *.  n ip.  ■ fn 


which  yields 


cov(*P 


n-1 


'n)-E{  r (n-m)At2(q^-E{q£))*Y(n-m)  £ At2(q*  -E(q?  , i ) > ^ \ s- . . it 2 Cq®- 

m-2  ii  m m,2  k.2  k i k i r.2  T 


n-1  m 


(6.1.13) 


-E(qjl)*  r (s-r)  £ At2(q*_1-E{q^_1)) 


n-1  s-1 

r r 


.4r.,_a  _a 


n-1  s-1  m r 


Z (n-m)(s-r)At'D(q“,q")+  £ £ £ £ (n-m) (s-r)At4D(q?  . ,q®  . )♦ 

nt-2  r«2  ™ T m2  r-2  k-2  n-7  k*’  t1-1 


r*2  k»2  p« 


n-1  s-1  r 


♦2  £ L E (n-«)(s-r)At4D(q“,q8 
*•2  r-2  p-2  ™ T-') 


p-2 


I 
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where  D(  ) represents  the  dispersion  of  the  enclosed  quantities.  But, 
as  we  have  already  discussed,  the  gradiometer  measurement  unit  is  com- 
pletely separated  from  the  accelerometers  and  gyros  units'  platform. 
Consequently,  it  is  reasonable  to  declare  that  the  gradiometer  and 
accelerometer  quantization  errors  are  not  correlated  and  therefore, 
the  last  teim  of  eq.(6.1.13)  is  dropped.  For  the  dispersion  of  the 
quantization  errors  an  exponentially  decreasing  correlation  function 
is  assumed  of  type (see  Denhard,  1977) 


(6.1.14)  D(ai,a.)u‘  e 


where  o the  variance  of  each  measurable,  t^-t^  the  sampling  interval 
and  t the  correlation  coefficient.  Therefore,  taking  into  account  that 
we  have  different  correlation  coefficient  between  accelerometer  unit(r 
and  gradiometer  unit(x  ),  then  the  variance -covariance  expressions  are 
given: 


a 


4 2 ^-1 

var(6Pn)»4t  c £ £ (n-m) (s- n)e 

“ m=2  n-2 


♦6tV 

g 


n-1 

r-1 

m*2 


5-  1 

Z (n-m) (s-n) • 
n«2 


m n 1 

L E e Tg 

k*2  p-2 


(6.1.15) 


. , n-1  s-1 

cov(6P  ,4P_)*At  a,  T.  T.  (n-m)(s-T)e 
n r 3 m-2  r-2 


VV 


♦«2 * 4c l 

a 


n-1 

L 

m*2 


s-1 

T.  (n-m)  (s-r)  • 
r-2 


m r 

T.  7 e 

k-2  r-2 


2 2 

where  o^,  the  accelerometer  and  gradiometer  variances  respectively. 

Eqs(6.1.15)  are  our  simulation  equations  to  find  the  order  of  magnitude 
of  the  quantization-induced  position  error. 
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6.2  Results 

Two  conputer  programmes  have  been  written  and  compiled  in  the  Siemens 
Computer,  the  first  to  examine  the  variance  case  and  the  second  the  co- 
variance  one.  These  two  programmes  are  listed  in  Appendix  A.  The  ground 
values  for  the  five  parameters  of  the  above  equations  have  been  chosen 
to  be: 

at  «0.1  sec. 

c2  -0.01  m2/sec4 

2 , c 
o “1 . E. 
g 

t -0.001  sec. 
a 

t -0.001  sec. 
g 

The  same  statistical  technique  used  in  simulation  I , is  followed  again 
for  the  quantization  error  simulation  studies.  Four  of  the  five  parame- 
ters have  been  kept  constant  and  at  the  same  time  the  fifth  one  assumed 
a range  of  possible  values.  The  computer  results  are  given  below: 


62 


1 


0.10  10'9  0.20  10'9 


0.30  10 


->■5 


.-5 


,-9 

-5 


0.40  10 


0.10  10  J 0.20  10  J 0.50  10  J 0.40  10 


-9 

,-5 


,-9 


-9 


0.50  10  0.80  10  0.11  10 
0.50  IO-5  0.80  IO-5  0.11  10 


.-4 


0.14  10'8  0.20  10"8 
0.14  10'4  0.20  10"4 


0.30  10 
0.30  10' 


Acceleration-induced  quantization 

error. The  first  rows  represent 

a =0.001  and  the  second  ones  o =0. 
a a 

(Covariance  case) 


0.50 

10'9 

0.60 

10'9 

0.70 

IO'9 

0.80 

IO'9 

0.50 

10'5 

0.60 

IO'5 

0.70 

10'S 

0.80 

IO'5 

0.14 

io'8 

0.17 

10‘8 

0.20 

IO*8 

0.23 

IO'8 

0.14 

10'4 

0.17 

10'4 

0.20 

io'4 

0.23 

io-4 

0.26 

10'8 

0.32 

10'8 

0.38 

,0'8 

0.44 

io'8 

0.26 

10'4 

0.32 

10'4 

0.38 

io'4 

0.44 

io-4 

0.40 

10-8 

0.50 

IO-8 

0.60 

IO'8 

0.70 

10‘8 

0.40 

1 

10'4 

0.50 

io'4 

0.60 

io'4 

0. 70 

io'4 

o.ss 

10'8 

0.70 

,o-8 

0.85 

IO'8 

0.10 

10-/ 

0.55 

IO-4 

0.70 

io-4 

0.85 

io'4 

0.10 

IO'3 

0.91 

io'8 

0.11 

1C'7 

0.13 

IO"7 

1 

0.91 

io-4 

0.11 

io'3 

0.13 

IO'3 

c4 

f 0.14 

.7 

1C 

0.17 

_7 

’0 

0.14 

IO'-5 

0.17 

IO"3 

! 0.20  :0’7 
i 0.20  ’O'3 


1 117,0 


At-a)O.Qlsec  jbjO.isec  , |c)0.5sec  M)1sec 


ja)0.1sec  b)1sec  cJ5sec  dllOsec 

0.6400 

io'10 

0.6400 

IO'6 

C.4000 

IO'3  |o.6400 

10'“ 

i 0.2 

( 

2 

10 

20 

0.3240 

IO"9 

0.3240 

IO'5 

0.2025 

IO'2  jO. 3240 

10'' 

j 0.3 

3 

IS 

30 

0.'340 

IO"9 

0.7840 

io"6 

0.4900 

IO'2  p. 7840 

io'1 

0.5 

5 

25 

50 

0.2304 

io'8 

0.2304 

io'4 

0.1440 

10' 1 10.2304 

j 

1 

10 

50 

100 

0.9604 

IO'8 

0.9604 

io'4 

0.6003 

IO'1  jo.  9604 

; 

20 

100 

200 

0.3920 

io'7 

0.3920 

io'5 

0.2450 

|3.9200 

3 

30 

ISO 

300 

0.888C 

IO-7 

0.8880 

io"3 

0.5550 

iS  .8800 
l 

1 4 

40 

200 

400 

0.1584 

IO'6 

0.1584 

io"2 

0.9900 

jl  5 .840 

5 

50 

250 

500 

0.2480 

10'6 

C.2480 

IO"2 

1 .S500 

124.800 

4 

6 

60 

300 

600 

0.3576 

IO'6 

0.3576 

IO'2. 

22.3500 

[35.7  60 

0.10  io'7  0.20  IO"7 

0.30  IO'7 

0.40  IO'7 

1 O.SO  IO"7 

0.80  IO'7 

0.11  IO'6 

. 1 

| 0. 14  IO'6 

0.20  IO'6 

I 0.30  IO'6 

0.50 

IO'7 

0.60 

IO'7 

0.70 

.7 

10 

0.80 

10 

0.14 

IO'6 

0.17 

IO'7 

0.20 

.7 

10 

0.23 

10 

0.26 

IO'6 

0.32 

IO'6 

0.38 

to'6 

0.44 

10 

0.40 

IO'6 

0.50 

IO'6 

0.60 

10'6 

0.70 

10 

1 0.55 

IO'6 

0.7o 

IO'6 

0.85 

IO'6 

0.10 

10 

0.91 

IO'6 

0.11 

IO'5 

0.13 

10 

0.14 

io'5 

0.17 

10 

0.20 

10' 

Gravity  gradients-induced  quantization  error (the  same  output  matrix 
from  IE  -1000E) . The  same  results  for  t. (Covariance  case) 


63 


o.to  10'7 

0.20  10' 7 

0.30  10'7 

0.40 

10'7 

0.50 

IQ'7 

0.60 

.o'7 

0.70 

to'7 

0.80 

ID'7 

0.10  10'7 

0.29  10'7 

0.S6  10'7 

0.91 

10-7 

0.13 

10'6 

0.18 

10'6 

0.23 

10"6 

0.29 

10'6 

0.50  10'7 

0.80  10'7 

0.11 

10'6 

0.14 

10'6 

0.17 

10'6 

0.20 

10-6 

0.23 

10-6 

0.86  to'7 

0.1?  10'6 

0.27 

JO'6 

0.40 

10'6 

0.55 

10'6 

0.71 

1C"6 

0.89 

10'6 

0.14  10'6 

0.20 

10'6 

0.26 

10'6 

0.32 

,0'6 

0.38 

to"6 

0.44 

10*6 

i 

0.33  10-6 

0.S5 

to'6 

0.80 

10'6 

0.11 

10'5 

0.14 

10'5 

O.18 

10-5 

| 

I 0.30 

10'6 

0.40 

10'6 

0.50 

10'6 

0.6C 

10"6 

O.'O 

1c'6 

J 

0.90 

10'6 

0.13 

10'3 

0.18 

10'5 

0.24 

>o'3 

0.30 

-O'5 

0.55 

10'6 

0.70 

10'° 

C.65 

10'5 

O.iC 

_c 

1C  " 

0.20 

10*3 

0.27 

10'5 

0.3t 

!0'5 

0.45 

.o's 

Correlation  coefficient- 

■induced  quan- 

! 0.91 

10'° 

0.11 

10'5 

0.13 

10' 3 

tization  error. 

The  first  rows  repre- 

1 

0.38 

10'S 

0.50 

1c'3 

0.63 

io's 

sent  t =0.002  and  the  second  ones  j 

, a . 10.14  lc'"'  0.1:  1C'6 

t =1  sec.  (Covariance  case) 

a J 0.65  10  ' 0.6.  10  3 

i — — ■ 

! 

' 0.20  10  5 

’■0.11  10'4 
! 


6.3  Discussion 

For  the  given  results  the  following  comnents  summarize  their  meanirg: 

a)  comparing  the  errors  commited  by  the  gradiometer-aided  navigation 
system  with  those  coining  from  the  quantization  studies,  it  can  be  cle- 
arly seen  that  the  latter  errors  constitute  a very  small  Quantity . The 
results  can  be  generally  considered  as  satisfactory  in  view  of  the  fact 
that  the  approximation  formula  used  for  the  inertial  acceleration  compo- 
nents approximates  grossly  the  reality  using  only  three  points.  If  moie 
terms  in  the  Stirling’s  formula  are  taken  into  account,  then  tit  quanti- 
zation error  studies’  results  would  be  effectively  reduced. 

b) the  sampling  interval  At  is  the  worst  contributing  error  factor  in 
the  quantization  error.  Decreasing  the  operation  of  the  system  to  the 
order  of  10(ten),  then  the  quantization  errors  increase  up  to  the  10'\ 

c) for  the  rest  four  parameters  of  the  system  only  the  accelerometer  va- 
riance causes  changes  in  the  quantization  error  budget.  For  exajple, 
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decreasing  the  variance  of  the  measured  acceleration  10  times,  then 

2 

the  quantization  errors  increase  141  to  the  10  . 

d) when  the  parameters  og,  x&,  xg  undergo  their  range  of  changes,  the  qu- 
antization errors  remain  unaffected.  For  that  reason  we  have  listed 
only  the  a contribution  with  the  understanding  that  the  rest  two  pa- 
rameters  give  identical  results 

e)  comments  b)  and  c)  are  also  valid  for  the  covariance  case  which  is 
included  herein  for  instructive  purposes. 

In  the  next  pages  we  give  some  representative  nomographs  to  picture 
briefly  the  quantization  error  studies  results. 


Atlsec) 
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7.  Simulation  II:  The  General  Error  Model 
7.1  Accelerometer  Error  Studies 

•It  is  well-known  that  the  accelerometer  frame  is  materialized  by  the 
three  input  axes  of  the  on-board  accelerometers.  Since  it  is  instrumen- 
tally  impossible  to  direct  three  axes  so  as  to  construct  an  orthogonal 
frame,  the  accelerometer  frame  is  finally  a non-orthogonal  or  quasi- 
orthogonal  frame.  Consequently,  the  fact  of  measuring  the  apparent 
acceleration  components  along  a non- orthogonal  frame  should  be  seriou- 
sly taken  into  account . Having  corrected  the  sensed  acceleration  for 
accelerometer  non-orthogonality,  then  it  refers  to  the  actual  platform 
frame.  A transformation  which  takes  the  acceleration  signal  from  the 
actual  platform  frame  to  the  ideal  one,  is  the  next  step  to  be  accom- 
plished. All  gyro  misfunctions  are  included  in  the  aforementioned  trans- 
formation. As  soon  as  the  apparent  acceleration  refers  to  the  ideal 
platform  frame,  which  in  our  case  coincides  with  the  navigation  frame, 
then  its  components  can  enter  the  general  equation  of  inertial  naviga- 
tion. 


Taking  into  consideration  what  is  discussed  above,  the  apparent  acce- 
leration signal  transformation  could  be  illustrated  by  the  general  re- 
presentation 


(7.1.1) 


c* 


1 Pa  a 
C Aa 
a a 


where  A^  , Aa  represent  the  acceleration  signal  coordinatized  in  the  na- 
vigation and  accelerometer  frames  respectively  and  C represents  the  di- 
rection cosine  matrix (from  where  the  notation  comes)  which  transforms 
the  frame  indicated  by  the  subscript  to  that  indicated  by  the  super- 
script. Pa  and  P1  denote  the  actual  and  ideal  platform  frames  respe- 
ctively. 


Since  we  are  using  very  often  transformations  of  type  affine  and  simila- 
rity ones,  it  would  be  helpful  to  define  them  from  the  beginning 
a) the  group  of  the  affine  transformations  can  be  represented  by  a ro- 
tation matrix  H plus  a transformation  vector  t,that  is 


(7.1.2) 


T(u)"Hu*t 
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(see  Grafarend  and  Schaffrin,  1976) . The  affine  transformation  preserves 
Euclidean  parallelism,  straight  lines  are  transformed  into  straight  li- 
nes and  planes  into  planes. 
b)when  the  two  relations 

H - XR 

(7.1.3) 

R_1R  -I 

hold,  then  the  affine  transformation  group  is  called  similarity  transfor- 
mation group  and  under  that  ratios  of  distances  and  angles  are  preserved. 

In  our  analysis,  skew-svmr.etric  matrices,  denoted  by  R^,  are  very  often 
used  to  transform  two  misaligned  orthogonal  coordinate  frames  into  eacr. 
other  in  case  the  misalignment  angles  are  considered  small.  It  helps  in 
the  understanding  of  what  follows  to  note  that  skew -symmetric  matrices 
are  always  transformed  under  the  similarity  group. 

Now,  we  shall  try  to  determine  the  two  transformation  matrices  involved 
in  eq. (7.1.1)  taking  into  account  the  error  sources  which  cause  them 
to  depart  from  the  identity  matrix. 

Pa 

1,  C - transformation 
a 


The  transformation  between  the  quasi  -orthogonal  accelerometer  frame  a and 
the  orthogonal  actual  platform  frame  Pa  is  a "small  angle"  transformation 
parametarized  by  the  small  angle  rotations  connecting  the  two  frames.  This 
transformation  is  treated  in  many  textbooks  in  detail  and  it  will  not  be 
further  considered  herein(for  a discussion  see  Britting,  1971,  p.39).  Ta- 
king into  account  the  angles  definition  depicted  in  Fig.  (9),  we  write 
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Fig.  (9):  Actual  platform,  accelerometer  frames  geometry. 

The  only  problem  on  which  we  like  to  draw  attention  concerning  the  above 
transformation  is  that  the  angles  inside  the  matrix  are  very  small  as 
well  as  the  off  diagonal  terms  are  not  equal,  e.g.e^.,  f The  justi- 
fication for  that  comes  from  the  non -orthogonality  of  the  accelerometer 
frame.  Consequently,  a small  rotation  about,  say,  the  a^-axis  first  and 
then  about  a^  and  a„  will  not  be  sufficient  to  bring  the  accelerometer 
frame  in  coincidence  with  the  actual  platfe*,.:.  We  note,  finally,  that  the 
six  angles  depicted  in  eq.(7.1.4)  can  be  measured  by  well-known  alignment 
techniques . 


C „ - transformation 
Pa 


As  we  said  before,  the  C -transformation  is  by  far  the  most  critical 

Pa 

operational  procedure  in  the  whole  navigation  systems  analysis.  Deeply 
thinking,  what  is  written  in  the  literature  known  to  the  writer  could 
be  considered  as  a mess  as  far  as  thus  transformation  is  concerned.  Con- 
sequently, we  feel  that  it  is  our  turn  to  put  things  into  an  order  by 
making  from  the  very  beginning  the  following  statements: 
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a)  When  the  moving  vehicle  is  at  the  starting  point  we  have  to  decide 
which  will  be  the  computation-navigation-ideal  platform  frame.  In 
our  analysis,  the  decision  was  taken  in  favour  of  the  Greenwich 
orthogonal  frame  and  it  was  inscribed  on  the  moving  platform  ever- 
af ter . 

b)  In  order  to  have  at  any  time  instant  the  ideal  platform  frame  paral- 
lel to  the  navigation  frame,  the  former  is  commanded  to  the  earth’s 
rotation. 

c)  The  following  statement  has  no  impact  on  the  mathematical  analysis 
of  the  problem  under  consideration,  but  it  has  to  be  made  in  order 
to  give  rigor  and  clarity  to  the  general  concept:  an  inertial  coor- 
dinate frame  is  somehow  and  somewhere  inscribed  on  the  moving  plat- 
form and  we  refer  to  it  when  we  postulate  that  the  platform  rotates. 

The  first  idea  to  be  accomplished  is  to  materialize  such  an  inertial 
frame  by  a set  of  three  single-degree-of-freedom  comoving  but  inertially 
stabilized  gyros.  Otherwise,  who  can  insist  on  saying  that  the  mo- 
ving platform  is  inertially  rotating? 

d)  The  command  for  platform  rotation  equals  to  the  earth’s  rotation 

is  injected  to  the  gyros  which  drive  then  the  platform  accordingly. 

But  since  the  gyros,  like  all  other  instrumental  units,  are  burdened 
wj.th  a variety  of  serious  errors  e.g.  gyro  drift,  they  have  to  be 
plugged  into  the  C.  - transformation. 

e)  At  the  starting  pofnt,  the  platform  frame  is  set  to  be  parallel  to 
the  navigation  frame.  Of  course,  this  is  by  no  means  true  and  thus 
an  initial  misalignment  is  everafter  present. 

Let  us  now  proceed  in  determining  the  discussed  transformation.  As  it 
is  shown  in  Fig.  (10),  at  the  starting  point (t=0)  the  actual  platform 
frame  Pa  has  a small  initial  misalignment  with  respect  to  the  ideal 
platform  frame  P1  due  to  the  errors  in  the  alignment  procedure.  These 
two  frames  seen  at  another  time  instant  t have  already  changed  their 
respective  attitude  due  only  to  the  inabilities  of  the  gyTOs.  The  actu- 
al platform  frame,  besides  its  initial  misalignment,  nas  already  got 
another  small  angle  distortion  , time  dependent  one,  denoted  by  the  ^ 
angles  6^ (i=X,Y,Z) . At  any  time  instant  t,  the  actual  platform  frame 
can  be  linked  to  the  ideal  platform  frame (or  for  that  matter  to  the 
navigation  frame)  with  the  general  transformation: 
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(7.1.5)  P1  " ^ “ V V* 

Pt  lo  1 
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where  the  R - transformation  relates  the  actual  platform  frame  at  time 
instant  t ^ to  that  actual  platform  frame  at  the  starting  point  and 

pi 

the  Rr  - transformation  is  the  initial  misalignment  transformation. 


Fig.  (10):  Platform  frames  variations. 

pi 

The  R - transformation  is  a constant  matrix  transformation  and  pre- 
sents  tQ  no  dificulties.  According  to  the  given  angles  definition,  we 
could  write(see  Fig. (10)): 


1 -IDj  my 

pi 

(7.1.6)  Ft  » m,  1 -mv 

p*  *■  A 

to 

“X  1 


Let  us  now  abstractly  denote  the  R 0 - transformation  by 


1 -6. 


(7.1.7)  R ® - fi„  1 

pa  L A 
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where  the  elements  of  the  matrix  represent  small  angles  or  rotation? 
due  to  gyro  errors.  Let  us  therefore  analyse  them: 
a) the  gyro  frame  is  constructed  by  the  three  spin  axes  of  three  on- 
board mounted  gyros.  They,  generally,  instrument  a non-orthogonai  or 
qua si -orthogonal  frame.  Consequently,  the  transformation  between  the 
quasi -orthogonal  gyro  frame  and  the  actual  platform  frame  reads. 


(7.1.8) 


where  the  $’s  represent  small  misalignment  angles  and  the  same  com- 
ments as  in  the  accelerometer  case  apply  to  the  gyro  non-orthogonality. 
b)nov,  we  demand  from  the  gyros  to  command  the  platform!  with  the  earth’s 
rotation,  but  since  the  three  gyros  have  generally  different  scale 
factor  uncertainty,  then  the  signal  for  the  respective  rotation  is  laj- 
sified.  The  gyro  scale  factor  uncertainty  matrix  is  expressed  as 

Ux  ° ° 


(7.1.9)  U - 0 Uy  0 
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where  U^.Uy  and  U-,  represent  the  X,Y  and  Z gyro  scale  factor  uncertainty 
respectively 

c)taking  into  account  eqs(7.1.7),  (7.1.8)  and(7.1.9),  we  find  the  error 
angle  6 to  be  expressed  as: 


(7.1.10) 


where  w^,Wy,u)z  the  earth’s  rotation  components.  From  the  last  equation, 
it  is  clearly  seen  that  the  gyros  channel  a signal  for  the  earth’s  ro- 
tation to  the  actual  platform  frame  plus  an  error  rotation  due  to  the 
gyro  non-orthogonality  and  scale  factor  uncertainty.  This  error  signal 

is  proportional  to  the  applied  rotation,  in  our  case  the  earth’s  rotation. 
Pa 

Consequently,  the  R to-transformation  can  be  analvticallv  written  as 
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The  off-diagonal  terms  of  the  above  matrix  are  time  dependent  quantities 
and  describe  that  as  the  gyro  frame  drifts  changing  its  angles  of  non- 
orthogonality, then  the  attitude  of  the  actual  platform  frame  is  affe- 
cted. In  case  in  which  the  gyro  frame  is  orthogonal,  it  is  not  drifting, 
has  no  scale  factor  uncertainty  and  in  absence  of  initial  misalignment, 
then  the  actual  platform  frame  is  nothing  else  but  the  ideal  platform 
frame.  It  is  therefore  seen  that  the  time  increasing  gyro  drift  causes 
the  above  off-diagonal  terms (those  inside  the  parentheses)  to  exist. 

Taking  into  account  eqs  (7.1.5),  (7.1,6)  and  (7.1.11),  we  find  the  ge- 
neral transformation  taking  the  accelerometer  signal  from  the  misaligned, 
non-orthogonal  accelerometer  frame  to  the  navigation  earth-linked  frame 
to  read: 
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Now,  since  the  apparent  acceleration  components  are  measured  by  three 
actual  accelerometers,  it  is  logical  to  assvme  an  overall  accelerome- 
ter error  model.  Taking  into  account  that  each  accelerometer  has 
its  own  scale  factor  uncertainty,  bias  and  random  uncertainty,  a gene- 
* ral  accelerometer  error  model  could  be  expressed  as  (Britting, 1971 , 

Denhard,  1977): 
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where  (i=X,Y,Z)  indicates  the  apparent  acceleration  output  signal 
of  the  quasi -orthogonal  accelerometer  frame 

AjJ1(i=X,Y,Z)  indicates  the  apparent  acceleration  as  an  input  in 
the  accelerometer  frame 

a^  (i=X,Y,Z)  the  accelerometer  scale  factor  uncertainty 
b|  (i=X,Y,Z)  the  accelerometer  bias  and 
u|  (i=X,Y,Z)  the  accelerometer  random  uncertainty. 

As  it  is  easily  seen,  eq.(7.1.13)  is  an  affine  transformation.  Needless 
to  say  that  the  left  hand  side  of  this  equation  is  the  acceleration 
signal  to  be  transformed  to  the  navigation  frame,  as  per  eq .(7.1.12), 
in  order  to  be  used  for  the  simulation  studies. 

7. Z Gradiometer  error  studies 

In  our  inertial  navigation  platform,  the  spherical  gradiometer  employed 
has  been  developed  and  tested  in  M.I.T.  Each  instrument  has  the  capabi- 
lity of  measuring  two  independent  gravity  gradients  and  therefore,  three 
of  them  could  furnish  the  whole  gravity  gradient  tensor  plus  a redundant 
gradient  indicating  accuracy.  As  it  is  intuitively  understood,  the  mea- 
suring process  is  quite  complicated  due  to  the  inherent  electronics, 
but  in  principle  the  following  fundamental  ideas  are  very  helpful: 
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a) Going  into  the  very  beginning  of  the  gradiometer  unit,  what  is  really 
sensed  and  measured  is  nothing  else  but  rotations  of  the  float  with  re- 
spect to  the  stabilised  housing.  These  rotations  are  sensed  by  a set  of 
electronic  axes  and  applied  back  to  restore  the  initial  float  attitude. 
The  measurements  of  these  rotations  represent  measurements  of  gravity 
gradients.  Taking  into  account  the  most  general  case  in  which  the  elec- 
tronic frame  is  non-orthogonal  and  is  slightly  misaligned  with  respect 
to  the  float  frame,  then  these  effects  have  an  error  influence  on  the 
measurements  of  the  gravity  gradients  which  must  be  anyway  compensated. 

b) The  float  frame  has  a certain  prescribed  orientation  with  respect  to 
the  axes  of  principal  moments  of  inertia.  But  due  to  various  reasons, 
e.g.  inability  in  locating  for  perfect  the  axes  of  principal  moments 
of  inertia,  the  float  frame  is  thus  considered  to  be  slightly  misalig- 
ned with  respect  to  the  ideal  float  frame. 

c) As  all  instrumental  packages  so  the  gradiometer  one  has  its  own  instru- 
mental axes  along  of  which  the  gravity  gradients  are  measured.  Our  grc- 
diometer  package  frame  is  the  ideal  float  frame  into  which  gravity  gra- 
dients measured  by  a non-orthogonal  electronic  frame  must  be  finally 
transformed. 

It  is  now  clear  that  on  each  float  four  different  sets  of  coordinate 
frames  exist.  These  frames  are: 

1.  Principal  moments  of  inertia  coordinate  frame (P^) 

2.  Gradiometer  measurement  unit  coordinate  frame (G^) 

3.  Actual  float  frame (F^)  and 

4.  Electronic  frame (E^) 

We  have  to  remark  that  the  P.  and  G.  frames  are  invariant  from  tne  gr3- 

li 

diometer  configuration,  but  the  frames  and  E^  do  depend  on  that  as 
per  Fig. (11). 

In  order  to  get  the  expressions  of  the  grt.  Tty  grauients  referred  v.o 
the  G^-frame,  we  proceed  as  follows: 

a)  from  the  general  gradiometer  torques  equations  find  their  respective 
ones  expressed  in  the  G^,  F^,  E^-frames,  for  the  torque  equations  are 
given  in  the  P^-frame. 

b) find  the  relations  between  torques  in  the  electronic  frame  and  gra- 
vity gradients  expressed  in  the  gradiometer  measurement  unit  frame  as 
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well  as  instrumental  non -orthogonalities  and  misalignment  angles. 

c) find  the  expressions  for  the  gravity  gradients  in  the  -frame  with  re- 
spect to  the  torques  in  the  s-une  frame. 

d) get  the  final  relationship  between  gravity  gradients  in  the  electronic 
frame  and  those  referred  to  the  gradiometer  measurement  unit  frame. 


As  soon  as  the  steps  a)-d)  have  been  carried  out  and  the  gravity  gra- 
dients in  the  -frame  have  been  got,  then  we  have  to  consider  a gene- 
ral miscrientatior.  of  that  frame  with  respect  to  the  inertial  one  to 
which  the  gravity  gradients  finally  refer.  Then,  the  gradient  tensor 
after  this  very  lengthy  procedure  could  be  used  in  the  final  equation 
of  inertial  navigation  if  anc  only  if  a special  coordinate  transforma- 
tion is  applied  to  "switch”  the  gradients  from  the  inertial  frame  into 
the  operational  earth- fixed  navigation  frame. 


Strictly  speaking,  the  mentioned  navigation  frame  cannot  be  used  since 
the  coordinate  differences  referred  to  it  cannot  be  integrated.  The  rea- 
son for  that  is  that  such  a frame  is  affected  by  time-like  misclosures 
due  to  polar  motion.  These  misclosures  have  been  computed,  but  in  terre- 
strial navigation  applications  are  to  be  safely  neglected (Doukakis, 1978) . 


As  it  is  seer,  from  Fig. (11 J,  the  moments  or  torques  measured  along  P1  and 
P?  principal  moments  of  inertia  axes  are  given (Trageser, 1975) : 


— ^lgp  P 

‘ 2 j 

(7.2.1) 

M2 

41gp.p, 

1 J 

AI 

*Ip2P2~Ip3P 

Eqs(7. 

2.1) 

will  now 

axes. 

Let 

us  first  ■ 

the  (P2P j) -plane (the  same  analysis  is  applied  in  all  configurations  by  a 
sinple  permutation  on  the  indices).  The  axes  ore  depicted  in  Fig. (12). 
From  elementary  plane  vector  calculus  we  get 


(?.2.2)  *^2cos45','^3sln4^ 

2y  *-?2sin4S*PjCos4^ 

Representing  the  gravity  gradient  field  by  r,  we  could  write 


(7.2.3)  r 


'P2Vp2P2*?3?3*P3P3 


+?2?3*P2P3+?!3P2«P,P 


3 2 
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Fig.  [12):  Axes  configuration  in  the  (P7P3)-plane  of  the  XY-float. 

where  the  g’s  represent  gravity  gradients,  Now,  each  gradient  can  be 
written  as 

(7.2.4)  g..  -t.(r.J) 

and  thus  we  get 

G G P P P 

"T'^YY'8XX^  >m  "T_(:8PjPj'8P2P7)cos2‘4S  'gP5P2sin2'45 

■ ^7-2-5)  =-Sp3P2 

where  the  letter  over  the  gravity  gradients  indicates  the  frame  which  they 
refer  to. 

In  view  of  eqs(7.2.1)  and  (7.2.5)  as  well  as  the  XY-float  configuration 
schematic,  we  get 

Mz  * T “ f*yy'gxx-) 

where  it  is  understood  that  the  gravity  gradients  refer  to  the  same  coor- 
dinate frame  as  the  measured  moments,  then  the  gravity  gradient  super- 
script is  dropped  for  simplicity.  If  eq(7.2.4)  is  applied  to  the  other 
planes,  then  we  can  get: 


‘ ^7(82X*8ZYJ 


81 


G . 7 

(7.2.6)  My  — ^(gzx*g2Y)  XY-float 

“z  ‘ ^Y^SJCX3 

By  a simple  permutation  of  the  indices,  we  obtain  for  the  YZ  and  ZX- 
float  configurations  the  following: 

\ m ^(«22-8yY) 

G AT 

(7.^.7)  My  = ^7(gxy+gxz)  Y2-flcat 

k > 

k a^7(Sv;+8«) 

G . * 

My  ''^7^Sy;'*'SyX1  ZX-£loat 

G A ’ 

M-_  " -7(8y2+8yx) 

£oas  (7.2.7)  express  the moments  coordinate zed  in  the  gradiometer  measure- 
ment unit  frame  and  hold  as  they  stand  for  the  actual  float  frame  as 
well  as  the  electronic  frame  changing  only  the  superscript  G by  F and 
£ respectively. 


We  have  now  to  find  the  relation  between  the  torques  measured  in  the 
electronic  frame,  being  the  gradiometer  sensor  for  gravity  gradients, 
and  the  gravity  gradients  coordinatized  in  the  gradiometer  measurement 
unit  frame  taking  into  account  instrumental  non-orthogonality  and  mis- 
alignment. First,  let  us  consider  >.he  case  in  which  the  electronic  fra- 
me is  a nor. -ortho gone!  one  and  also  misaligned  with  respect  to  the  actual 
float  frame(see  Fig. (13)).  As  we  have  already  explained  in  the  accelero- 
meter studies,  the  matrir  which  transforms  the  moments  of  the  electronic 
fram  to  those  of  the  float  frame  is  exactly  as  per  eq.(".1.4).  Therefore 
we  write 

F E E E 

MX  * "x^YzV^Z 

(7.2.8)  fly  -ex2  k -V*zA 

c E E E 

'"z  -eXYVeYxVM2 
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Fig.  (13):  Relation  between  actual  float  and  electronic  frames. 


Solution  of  eqs (7 . 2.8)  with  respect  to  the  electronic  torques,  disregar- 
ding products  of  small  angles,  gives 
E F F F 

“V0xzMremMz 

(7.2.9)  ^ =-eY;VV9YX^; 

F.  r F F 

M:  m*zr\  -9:x'V  *M: 


and  since  the  electronic  torques  are  equal  in  magnitude  and  opposite  in 
sign  to  the  torques  caused  by  the  gravity  field,  we  get 
E F F F 

*x  *-VexzV6xYq: 

(7.2.10)  ^ -e^-Ry-eYxR. 

E F F F 

fIz  "ezA-ezA-*z 

where  the  bar  over  the  torques  indicates  that  these  torques  come  from  the 
physical  entity  of  gravity. 

Now,  let  us  transform  the  gravity  gradients  which  refer  to  the  gradiome- 
ter  measurement  unit  frame  into  those  referred  to  the  actual  float  frame. 
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Taking  into  account  that  the  float  frame  is  misaligned  with 
misalignment  angle  about  the  G^-axis 


then,  the  gravity  gradients  sensed  by  the  XY-float  configuration  reac. 
F G G G 

gxx ' *xx‘2*xyV28xy*3 

F G G G 

*yy  ’ &yy*28yzV28xyv3 

C7-2-11)  F G G G G G 

8zx  * Szx^V^xx-Szz^zy^ 

F G G G G G 
gZY  * SZY*  ^8ZZ~8YY^  ’ 1 *8XY*2  -*ZX*3 


V>e  rewrite  now  tne  iirst  three  equations  from  eqs(7.2.<)  in  t.ie  ac^ua, 

float  frame  as 

F l\  F F 
H " T^zx^zy' 

- ai,f  r 


(7.2. 12a)  My* Z 


F F F 

M,  1 %t*YY*gXx'J 


Solution  of  the  above  equations  with  respect  to  the  gravity  gradients 
gives 


F f 


M 


8ZX*8ZY*  ^ \ 

(7.Z. 12b)  gzx4;Y-TFV 
F F ..  F 
gYY'8XX  ~7  S1Z 

Combining  eqs (7.2.10),  17. 2.11)  and, 7. 2. 12),  we  sec  that  the  X-channel 
electronic  torque  assumes  the  form 


E.jGG  GG  G GGG  (■  G GG 

ftX  ’ ^^^ZX-^1+(Sxx-g;;)iVgZY^+gZY+Cgzz-R-A.H1*gXYe;-S-;^J-^(g:x-gVA. 

GG  GGG  GG  GGG  G GG 

♦(gxx-gzz)'i'2*gZY*(gzz-gYY)<l'1^2-gZX’i,3-^(8YY-gXX*2SYZ*r28XY^-,8X^3+gZYv3J' 


J 


After  some  manipulation,  we  obtain  the  electronic  torques  for  all  three 
float  configurations  to  read: 


84 


Ew  at  G G GGG  GGG  GG  GG 

flX  “ T ('^ZX**ZY^+(gYY‘gZZ*gYX^1+(gZ2'gXX,'gW)l<'2*(gZX'*ZY)'('J*^ZX*e2r^8XZ* 


G G 

*tgYY'gXXJ  6xy} 

Ew  A]  G G GGG  GGG  GG  GG 

" T ((82X+gZY)*(8ZZ‘gYY_gYxHi  + C8xX‘8ZZ'gXy)l<'2+(gZY"gZX)’^+(8ZX"gZT)6YZ' 

G G 

* ^gYY~8XX^  SfX1 

E /y  G G G ^ ^ GG  GG 

NlZ  = T ^gXX'8YY)‘2gYZ'h'2gXZ'*G:+4gXYl^'(gZX+gZY5  9ZY_(gZX+8ZY)  ^ 


(7.2.13) 


,.YZ  il 


*X  “ T_ffg'l7'g:2)+4gY2^'2gZX|J2'2gYX'(3'(8XY+8XZ)eXZ'(gXY+8XZ)0XY) 


Ey-^CG  GG  GGG  GGG  GG 

«Y  ’ r(-  (S)CY+gXZ)  + (gXY’gXZ^r  ^gZZ_gXX+gZY^2+  ^XX^YY^YZ5^*  (gZZ_8YY)eY:* 


G G 

k(gXY+gXZ^6YX) 


Ey,  AT  G G GG  GGG  uuo 

»?  - ^UgXY+gx2)-(gx2-gXY)*1+^xx-gYZ'gzz)',,2+(gYY'gxx'gYZ:ih  'g:Z*gYYjeZY" 


GGG 


G G 

•(gXY+gXZ)6ZX) 


:T((£YZ*8YX)-(g;:-gYY-gzx)^+(gYx-gYz)'i,2+-gYY'gxx'gxz)'i3‘(gxx'g:z)fx:' 

G G 

+ (gYZ+gYX)9XY} 

E7y  at  G G G G G GG  GG 

‘T " ^^gZZ'gXX^2gZY*1*4g2X,,'2'2gXY,i'3-(gYX+gYZ)9YZ-(8XY+gY:)9YX) 

E7y  .tGG  GGG  GG  GGG  GG 

Mp  ^(-(gYX+gYZ^(gYY-gZZ+g2X)4,1  + (gYZ-gYX),,;2+(8XX-8YY+8XZ)',,3*(8YX+8YZ)e:Y+ 
G C 

*(gxx"gzz)8zx) 


G G GGG 


GGG 


Now,  eqs(7.2.7)  are  written  in  the  electronic  frame  as  follows: 


~7T  "x  * ^ZX+^ZY 

“SI  *'^zx"^:y 


2 ^Y  . E E 
TT  *Z  gYY  8 XX 


+ ‘ M, 
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Having  excluded  eqs(7.2.14b,f,ii , the  rest  six  equations  are  the.*1  solved 
in  view  of  the  Laplace  condition 

(7.2. is)  &xx*8YY*szz  — 4-„ko  ♦ zj' 

where  w the  earth’s  rotation  with  respect  to  the  inertial  space,k  the  uni- 
versal gravitational  constant  and  o the  density  of  the  medium  in  which 
the  navigation  takes  place.  The  solution  gives: 

E 4 t 2 2 2 .„5zx  J-.iZ, 

gxx  *'  u + TiT^V  1 


(7.2.16) 


!» - 

8X2*  ZT^2*^ 


2 


gyy  - ♦»£“) 

F F F E 

1 -fiXY  jJZ 

g'.’Z  7TiMX  _My  _MX  ' 


E 4,2  2 2 *d.X 

gz;  «-  jwkp-^  u,  - m(}..  -M*  ) 


If  the  electronic  torques  involved  in  the  above  equations  are  substitu- 
ted by  those  given  in  eqs(?.2.13),  then  a^ter  a tedious  manipulation  we 
get  the  relations  between  gravity  gradients  in  tne  electronic  franc  and 
gradiometer  unit  frame  to  read: 


4221^^^^^  G -C 

gxx-  7 ny°*  7 w + 7(«zz+«n'2*xx+68zx*2-6Sxy^--2eYX*exy*ex;+2ev:)8)n.-(eXI»e)C>  )gx 

-(6yX+2ey2)gX2) 
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i 

i 

i 


E.GGGGG  G G G 

*X^“T  (^2«X^2-2«YZ*2-2«XxV28YY^-tBW*eXZJ*XX*teXZ-eYX)*XZ^eYZ*eX»J^Z- 
G 

-(0xz*eYZ)gzz) 

E.GGGGG  G G G 

8XZ“  7 (‘2gXZ+i:8XY4'l  + 2gZZ'l'2*28X2'<'rZ8XX'J'2'(0XY+eX2)gXX*(6XY‘,'eyX)8XY*(9XY+eyi)8YY+ 

G G 

(7.2.17)  +(9XY+9XZ)gYZ*(8YZ'eXZ)gZZ) 

E 4 G G G 

gyY-  ? »*>♦  y « - j(-gxx-gzz*2gyY^g^r6gXY*3-(26XZ+9YZ*9YX+29XYJgXY-(9YZ*9YX)8YZ- 

-(26XZ*2exY)gxz) 

E^GGGGGG  G G G 

gYZ"  7('2gZY'f2gYY'vr28:z'*'l,'28XZ,i'l*28ZX4'3'2gZY4'3'(eYX'9XZ^8XX'^eYX','9X'YJgXl'*(6XI_6Y'XJ8XZ* 
G G G 

+ (9YZ*eXY)gYY+(9XZ'9XY^gYZ+(9XZ‘9YZ)8ZZ) 

E 4 . , . G G G G G G G 

gZZ"  J u - j(282Z'gXX'8YY'68^r6gXZ4'2+(0XZ"eYX‘eY:+eXY-gXY*(6XZ*i;(Y'gX:- 

G 

^6YX+8YZ)gYz) 


Eqs(7.2.17)  can  be  written  in  the  concise  matrix  form  of  eq.  (7.2. 18). 

Now,  the  G-frame  is  considered  to  be  misaligned  with  respect  to  the  iner- 
tial frame.  Denoting  by  superscript  I those  gravity  gradients  referred 
to  the  inertial  frame,  then  we  write 


(7.2.19) 


3 


1J 


. ik' 

k*11*l 


a-’-3ilgkl 


where  i,j=X,Y,Z  and  the  a’s  represent  the  element  of  the  transformation 
matrix  between  the  two  frames.  Analytically,  the  gradients  are  written: 


8xx“axxaxxExx+axxaXYgxY+axxaxzgx:*aXYaxxeYX+aXYaxY8YY+aXYax:8Y:+ax:“xxB:i'axzaXY8:Y* 

*axzax:gz2 

i 

8XY*aXXaYX8XX+aXXaYY8XY+aXXaYZgXZ+aXYaYX8YX+aXYaYY8YY*aXYaYZ8Y'Z*aXZaYX®Zi*d:\'YaY'.g7Y+ 

+axzavz8:z 

I 

8XZ°aXXaZX8XX+aXXaZY8XY+aXXaZZBXZ+a.XY-aZX8YXtaXYaZY8^*aXYaZ:EYZ"aX:a:X8a*JXZa:Y8:Y* 

+aXZaZZgZZ 


(7.2.20) 
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I 

*YY'^X^XX*^X^XY+^X^ZgXZ*aYYaYXgYX*aYYaYY«YY*aYyaY2gYZ+aYZaYX8ZX*aY2aYYgZY* 

j ♦aYZ*YZgZZ 

8Y2‘aYXa2X8XX+aYXaZYgXY<'aYXa228x2+aYYa2x8'rt*aYVaZY8YY+aYYa228Y2*8Y2a2X82X+aY2aZYgZY+ 

+aY2aZZgZZ' 

I 

8z2'aZXaZX8XX+aZXaZY8xY*aZXaZ28xZ+aZYaZXgYX+aZYaZY8YY+aZYaZz8YZ*aZZaZX8ZX+aZZaZYgZY+ 

*azzazzgzz 

where  all  gravity  gradients  at  the  right  hand  side  refer  to  the  G-frame 
and  the  superscript  G has  been  dropped  for  simplicity.  Since  the  misalign- 
ment between  the  G-frame  and  the  inertial  one  is  in  the  "small  angle" 
sense,  then  we  can  make  the  approximations: 

axx'1  • azz“^Y 

(7.2.21)  *xv’^Z  aYY*1  aZ Y_'CX 

aX2*’5  Y aYZ"5X  aZZ’1 

where  the  represent  small  rotation  angles  about  the  Gx,Gy  311(1 

G7  gradiometer  axes  respectively.  Combining  eqs(7.2.20)  and  (7.2.21), 


l c G G 

gXX“g  XX"  :MxZ*  z[-  Zg  XY 

I G G G G G 

gXY‘‘gXY*c,XgXZ‘tYgZY^:fgYY‘8XX') 

I G G G G G 

gXZ'gXZ'SXgXY<Y(gXX'gZZ)+SZgYZ 

I G G G 

8YY*gYY*2^XgZY*2^ZgXY 

I G G G G G 

gYZ’'gYZ*f:X(8ZZ'gYY)+f-Y8YX'SZgXZ 

1 G G G 

8ZZ”8ZZ‘2'xgYZ*2MxZ 


we  get: 


(7.2.22a) 


Eqs(7.2.22a)  will  now  be  written  in  matrix  form  for  further  reference  as: 
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Now,  the  gravity  gradients  refer  to  the  inertial  frame,  but  in  order  to  be 
used  as  those  entering  the  fundamental  equation  of  inertial  navigation, 
they  have  to  be  transformed  into  the  navigation  frame.  The  transformation 
matrix  between  the  selected  inertial  frame  and  the  earth-linked  navigation 
frame  is  parametarized  by  the  true  sidereal  time  9 and  the  coordinates 
of  the  polar  motion  x and  y referred  to  a specified  epoch.  This  transfor- 
mation is  given(Veis, 1962) : 


cosc 

-sir.3 


-xcoi.P-ysu.e 


sinf 

cos& 

-xsin6*vcos3 


-y  L1] 

A 


and  according  to  standard  literature  the  gravity  gradients  are  transformed 
as  follows: 
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where  the  index  over  the  gravity  gradients  again  indicates  the  reference 
frame  in  which  they  are  coordinatized.  The  manipulation  of  the  above  matrix 
equation  gives  for  the  six  gravity  gradients  of  interest  the  following: 

1 2 22222  2 

gXX"gXXC0S  e+gyYsln  6't'8ZZ*-x  cos  ®+y  sin  ®)~8xsin20-2gxzxcos  e-gx2ysin28*g7Yxsin26+ 

♦g72xysine+2gZY>,sin20 

I 2 2 1 

8XY*8XXcos0sin8','8XY^1"2sin  0)+8X2(Y-2ysin  9-xcos9sin0  -xcos^-gyySinecosO* 

2 2 2 2 2 
+gY-Cxsinecose-ysin29-xcos  9)+gzz(x  cos  9-xycos  9+xysin9cos8-y  sinScose) 

I 2 

8XZ=8XXxcose"8XY  Cxsine-*-ycos03  “8xz(x  cose_cos6+xysine)^8YYys;'-ne"8zz^xcose*sin0^  + 
(7.2.23a)  ♦g,y(xycos0+xysin9-sin6) 

1 ■>  2 2 2 2 2 
gYY=gXXsin  0+8XYs*n2e*gXZ^ys*n28'2xs:*n  8)+g^(2ycos  6-xsin0)',-gzz(x‘'sin  0+y  cos“8' 

-xysin9) 

1 ? 2 
gyZ*gxxxsin0*gXY(xcos6-ysin0)+gx7(sin8-x‘sin9+x>'cose)-gyYycose+gy2(xysine-y  cos0+cos6)* 

*8zz(yc°s8-xsin9) 

^2  2 
*zimW  -2gxYxy*28x:x+8YYy  -2gY:y+g:; 

where  all  gravity  gradients  in  the  right  hand  side  refer  to  the  navigation 
frame.  Eqs(7.2.23a)  written  in  matrix  form  read (see  next  page). 

As  we  have  already  discussed,  the  first  procedure  to  be  followed  in  trying 
to  simulate  the  fundamental  equation  of  inertial  navigation  is  to  express 
all  vectoi  quantities  to  the  same  coordinate  frame  and  particularly  to  the 
navigation  frame  in  order  the  results  to  be  referred  to  the  earth-linked 
frame.  As  far  as  gravity  gradients  are  concerned,  the  following  transforma- 
tions must  be  made  to  "channel”  the  signal  as  it  is  sensed  by  the  electronic 
float  axes  to  that  referred  to  the  navigation  frame: 
l)rewrite  eqs(7.2.18),  (7.2.22b)  and  (7.2.23b;  absractly  as 

(7.2.  t8)’  8y  * [A]  *[b] 

(7.2.22b)-  8ij  - [C]  8ij 
(7.2.23b)-  8ij  - [DJ  g.. 
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where  the  g’s  represent  the  column  matrices  of  the  gravity  gradients  and 
A,B,C,D  the  already  defined  transformation  matrices. 

2) derive  those  gravity  gradients  referred  to  the  navigation  frame  using 
the  following  formula: 

(7.2.24)  gAj  - [d]'1  gAj  -[A])) 

where  the  superscript  -1  indicates  the  inverse  matrix  operation.  The  last 
equation,  in  view  of  the  involved  matrices,  shows  clearly  that  the  sensed 
gravity  gradients  should  undergo  a very  tedious  and  lengthy  manipulation 
containing  non-orthogonality  effects,  instrumental  misalignment  etc.  in 
order  to  be  finally  used  in  the  simulated  navigation  equation  to  be  pre- 
sented next. 


7.3  The  simulated  navigation  equation 

In  view  of  eqs(2. 1 1) , (7. 1 . 12) , (7.1 .13)  and  (7.2.24)  the  complete  naviga- 
tion equation  which  can  estimate  the  instantaneous  geocentric  coordinates 
of  the  moving  object  with  respect  to  the  earth-linked  navigation  frame, 
is  written  in  matrix  form: 

(7.3.1)  - T(aAJn*b*u)-D'1C(B'1(g^-A))-2^RE-^xR-f^x((^xR) 

where  the  above  symbology  has  been  already  defined.  Now,  in  order  to  get 
out  of  the  navigation  system  some  indicated  numbers  or  order  of  magnitude 
of  errors,  certain  assumptions  have  to  be  made (of  course,  a rigorous  sta- 
tistical analysis  should  include  the  full  matrices  involved): 

a)  Since  the  navigation  system  is  simulated  for  its  performance  during  a 
very  limited  time  span,  or  as  a matter  of  fact  for  a few  seconds,  it  is 
reasonable  to  neglect  all  terms  containing  angular  velocity  or  acceleration 
of  the  earth  coordinatited  in  the  inertial  space.  Consequently,  the  last 
three  terms  in  eq.  (7.3.1)  are  for  our  simulation  studies  dropped  out. 

b)  The  matrix  T given  in  eq .(7.1.12)  contains  products  of  small  angles 
plus  some  single  terms.  Having  decided  to  keep  only  first  order  terms,  then 
the  matrix  T can  be  approximated  by: 
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(7.3.2) 


1 -(eX2+m2) 


m,+Ev-  1 

L \C 


-(my+e^)  mx+e2X 


eXZ+mY 

■(eyx’,iix5 

1 


higher  order  terms 


Combining  eqs(7.1.13)  and  (7.1.2),  we  get  for  the  apparent  acceleration 
components  to  be  used  inside  the  navigation  equation  the  following: 


a 

Ax 

a i 

i* 

1 

i 

i 

C .3.*J 

■S 

ii 

A:' 

i 

1 L 

f 1 ♦ax)A"Vo>;*Ux-  (cX2*m2)Ajn+ 


(m2*i_y-)Ax  * (1  +a^,)A^n+by+u,.-(eyX*m2lA2n  j + higher  order  tori.i 


- (mY*-c;.y)A^n+  ( VCZX^  Ay"’ *a  )A;r‘+b-’+u;  j 


The  above  equation  shows  clear!)’  that  tht  components  of  the  apparent  acce- 
leration in  each  channel  are  a mixture  of  all  three  encountered  accelera- 
tion components. As  it  can  be  also  seen,  each  channel’s  signal  is  merely 
composed  of  its  counterpart  acceleration  and  the  other  two  signals  multi- 
plied by  small  quantities (due  to  the  accelerometer  non-orthogonality)  are 
present  as  well. 

c)  Let  us  now  consider  the  case  of  the  gravity  gradients.  As  eq.{7.2) 
shows  the  determination  of  the  gravity  gradients  to  qe  used  in  the  navi- 
gation equation  requires  the  inversion  of  two  matrices,  namely  the  D 
arwi  S'1.  Since  we  have  first  to  derive  the  final  error  matrix  in  front 
of  the  gravity  gradients  referred  to  the  electronic  frame  ana  then  to 
make  the  assumptions,  the  mentioned  inversions  must  be  computed  by  hand. 
B-1  presents  no  difficulties,  but  D-1  possesses  certain  problems  ate  to 
its  elements’ complexity.  Owing  to  this  fact,  eq. (7.2.22b)'is  considered 
as 


Therefore 


gE  -A-BC-y  -A+BC*'DgN  «0gN 
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arel  finally 
(7.3.4)  gN  -n'V 

The  inversion  of  C_1  by  the  partitioning  method  gives  the  following 


matrix  approximated  only  up 

to  the 

first  order  terms: 

1 

2 S' 

0 

0 

0 

h 

k 

-h 

^Y 

0 

(7.3.5)  c’'  * 

S 

fx 

i 

0 

fY 

♦ higher  order  terms 

0 

2^z 

0 

i 

~zk 

0 

0 

h 

1 

_ r 

0 

0 

-2^ 

0 

2^ 

1 

i 

Therefore,  the  transformation (7. 3. 4)  yields  the  following: 


•N  c r.  s 4 2 

gXX  * gXX+gXY  ,26 

N £ £ E E . E S 4 2 

SXY  “ gXY+gXX  'e+gX2  -y+gYY  -e  + I(eXZ+9YZ)gZZ_(  I "k3  ' T “ } 

NEE  E 

gX"  * 8X7-gy2  * Ei"*" ('+<2~  Sy*  I(;3y7‘exzj)gzz 

(7.3.6) 

N E E g 4 2 

®yy  = gn'gXY  ,28'(  J "k2  ' 3 •“  ] 

N E . E 

gYZ  * gYZ+^  Sc~^1+  7(eXZ*0YZ)  JgZZ 

NEE  E 8 4 , 

gz:  “ g:z'gxz  '2x+gYZ  ,2y'(  T *kc  ' I 

As  we  can  see  from  the  above  expressions,  the  gravity  gradients  referred 
to  the  navigation  f:ame  are  equal  to  their  counterparts  sensed  and  mea- 
sured by  the  electronic  frame  plus  products  of  gravity  gradients  with 
small  quantities  e.g.  polar  motion  components.  .All  these  terms  in  view  of 

1) the  conclusions  drawn  in  simulation  I with  regard  to  gravity  gradients 

2)  the  small  magnitude  of  all  terms  but  the  single  gravity  gradients  and 

3)  the  objective  of  our  simulation  studies,  are  finally  neglected.  Con- 
sequently, the  gravity  gradients  which  enter  the  simulated  navigation 
equation  are  those  which  are  measured  by  the  electronic  float  frame. 


Recollecting  results,  the  equations  which  we  shall  simulate  assume  ts 
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following  forms: 


(7.3.7) 


X3  * at2(tUVA^"bX^-^^^(cXZ^)^^1^1(X2-X1)"CXv/Y2-Y1>4 

+GX2/Z2-Z1»2X2-X1 

Y3  ■ At2((mrcY2)^(UaY)4^VV(cYX^^,^/X2-X1)*G^.j(Y2-Y1)* 


*n«2-W*2-h 

h * At2(-(Vc^)AiNVc2X)^(1*az)^*b2^:*G2i*G2Xi  CX^X,)*^  (Y2-Y^)  + 


+Gz;|C-;-Z1))*222-Z1 


The  above  equations  contain  36  parameters,  namely: 

a)  time  span(1) 

b) apparent  acceleration  components (3) 

c) the  geocentric  coordinates  of  the  first  two  points (6) 

d) the  gravity  components (3) 

e) the  gravity  gradients (6) 

f) the  accelerometer  bias (3) 

g) the  accelerometer  random  uncertainty (3) 
f)the  accelerometer  non-orthogonality ^5) 

i) t'ne  initial  misalignment  angles  (5) 

j) the  accelerometer  scale  factor  uncertainty (3) 

The  same  procedure  followed  in  simulation  I will  be  repeated  to  see  the 
influence  of  these  36  parameters  on  the  derived  coordinates 

7.3  Results 


The  same  computer  prograime  used  for  simulation  I studies  is  employed 
to  accommodate  the  new  simulated  equations.  Only  the  results  for  the 
parameters  which  can  contribute  changes  into  the  system’s  error  budget  ^ 
are  listed  below. 


— — 
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var 

A 

O 

< 

(m2) 

4.0029 

0.0029 

0.0437 

5.0029 

0.0437 

5.6626 

4.0129 

0.0029 

0.0437 

5.0029 

0.0437 

5.6626 

4.0529 

0.0029 

0.0437 

5.0029 

0.0437 

5.6626 

4.1029 

0.0029 

0.0437 

5.0029 

0.0437 

5.6626 

4.2029 

0.0029 

0.0437 

5.0029 

0.0437 

5.6626 

4.5029 

! 

0.0029 

0.0437 

5.0029 

0.0437 

5.6626 

6.0029 

0.0029 

0.0437 

5.0029 

0.0437 

5.6626 

9.0029  0.0029 

5.0029 


var 

-tov*:  (m2) 

1.0000 

0.0000 

0.0002 

5.0001 

0.0002 

5.0034 

1.0400 

0.0000 

0.0002 

5.0001 

0.0002 

S.0034 

1.2000 

0.0000 

0.0002 

5.0001 

0.0002 

5.0034 

1.4000 

0.0000 

0.0002 

5.0001 

0.0002 

S.0034 

1.8000 

0.0000 

0.0002 

5.0001 

0.0002 

5.0034 

3.0000 

0.0000 

0.0002 

5.0001 

0.0002 

5.0034 

9.0000 

0.0000 

0.0002 

5.0001 

0.0002 

5.0034 

21.0000 

0.0000 

0.0002 

5.0001 

0.0002 

5.0034 

\ (m'/s 


vur.-cov. 


5.0000  0.0000  0.0002 

S.0001  0.0002 

5.0034 


var.-cov. : 

On2) 

S.0000 

0.0000 

0.0002 

5.0000 

0.0002 

5.0034 

S.0000 

0.0001 

0.0002 

5.0001 

0.0002 

5.0034 

5.0000 

0.0002 

0.0002 

5.0003 

0.0002 

5.0034 

5.0000 

0.0003 

0.0002 

5.0006 

0.0002 

5.0034 

5.0000 

0.0006 

0.0002 

5.0012 

0.0002 

S.0034 

5.0000 

0.0009 

0.0002 

5.0018 

0.0002 

5.0034 

S.0000 

0.0012 

0.0002 

5.0024 

0.0002 

5.0034 

S.0000  0.0000 

0.0002 

0.02-0.50 

5.0001 

0.0002 

5.0034 

5.0000  0.0000  0.0002 

5.0001  0.0002 

S.0033 


5.0000  0.0000  0.0003 

5.0001  0.0002 

5.0033 


5.0000  0.0000  0.0004 

5.0001  0.0002 

5.0036 


5.0000  0.0000  0.0005 

S.0001  0.0002 

5.0039 


5.0000  C.0000  0.0008 

5.0001  0.0002 

5.0045 


S.0000  0.0000  0.0014 

5.0001  0.0002 

5.0057 


7.4  Discussion 


Needless  to  note  that  the  results  of  the  new  simulation  are  almost  iden- 
tical to  those  derived  in  simulation  I and  therefore  the  same  comments 
are  also  applicable  here.  But  since  we  are  not  at  all  satisfied  with  the 
system’s  behaviour,  we  try  to  investigate  the  aided  navigation  system 
deeper  by  making  the  following  studies: 

a)  Investigation  on  the  contribution  of  the  sampling  interval  At: to  moti 
vate  our  discussion  we  remiai  again  that  the  inertial  acceleration  compo 
nents  have  been  approximated  by  the  Stirling’s  formula(see  eq. (5.2.35). 
In  approximating  derivatives  a dominant  error  source  is  the  input  errors 
themselves.  One  could  immediately  see  the  explanation  to  tout  looking 

at  the  mentioned  equation.  The  reciprocal  power  of  the  sampling  interval 
At  multiplies  the  true  values  as  well  as  their  errors  and  thus  the  algo- 
rithm magnifies  them  enormously,  "or  that  reason  we  investigate  the  cose 
in  which  At  decreases  in  order  to  see  how  space  traverses  of  such  ar. 
aided  navigation  system  behave.  The  results  are  listed  in  the  next  page. 

From  tr.e  given  results  it  is  obvious  that: 

1) The  navigation  system  becomes  a little  more  tolerable  as  At  decreases 
but  again  it  be  laves  badly 

2)  If  we  call  point  10  of  the  first  traverse  as  the  "break  point"  of  the 
navigation  system,  then  it  occurs  at  point  11  and  13  for  the  rest  t raver 
ses. 

Consequently,  by  making  At  very-  small  the  system  u es  not  bcliave  better 
except  some  very  smal I improvements . 

b)  Investigation  on  the  contribution  of  the  omitted  covaria:.  es.frcm  the 
giver  results  of  simulation  studies  I and  II  it  is  evident  that  the  sys- 
tem is  very  sensitive  to  the  variances  of  the  initial  coordinates  and 
especially  to  the  variances  with  vhich  the  geocen^.ir  coordinates  of  the 
second  point  are  known.  Therefore,  we  try  to  see  if  the  omission  of  >.hf 
respective  covariances  makes  any  changes  into  the  system’s  error  budget. 
The  formulation  of  the  new  investigation  follows  these  brief  lines: 
1)rewrite  the  general  simulated  equations  as: 


ISO* 
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*i  - At' ‘ ( A * CL.  ♦ (O,  (xi-rxi-2)*  Gxy  or,  ,-y,  -,)♦ 

1 *i-2,i-1  xi-2  XAi-2 , i-1  1 1 12  Mi-2,i-l  1 ' 1 - 

CjBj.2.,-,  CZi-l'2i-2))+2Xi-2‘Xi-1 

Y - fit2  ( A^.  ♦Gy  ♦ Gyx  (Xi_rXi-->+  Sy  (Yi-rYx-^* 

1 ’i-2,i-1  ^i-2  * i-2, i- 1 1 ‘ 1 “ 1Yi-:,i-1  1 1 1 - 

Y“i-2,i-1  * 1 1 - * " 1 1 


7j  “ Af  ( 


7 + G,  + C-x  (X ■•.i-xi«?)+  C-.Y  O'-  .-Y  ,)< 

-1-2, i-1  “i-2  “i-2, i-1  * ' 1 " 11  12 


“^i-J.i-l  1 w 1 “ 11 

2)call  y the  parameters  of  the  a!'C're  equation  except  the  coordinates  ana 
C the  desired  coordinates.  Then  for  the  first  point  one  gets: 

& ’s  !3t  ‘M-i 

tr  * a.:  at 

Cj  .■>  y\> 

since  var(X1,Y1,irX2,Y,,:,)-C.  Ec_  indicates  the  dispersion  matrix  of 
the  coordinates  at  point  5. 


point  4 


[S]“  A4>4+B4f5  * [A4  B4]Efc] 

|C- 

L 21 


2|'  - iA,  Bi.lr, 


C,  [_4  4J 


Vl‘  »4] 


Y , r 


y4]L“  Si 


SI 


X4 

0 I, 


M 


Sj 


[A4  Bf 


y* 


0 tC;i 
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point  5 


[S]  * WSVSS  * fs  BS  Ds] 
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£C  *Ps  B5  Ds] 


[As  Bs 


Vsl 

m 

\ 

0~ 

'c; 

A4  B4  y4 

C4 

_C3 

_°  1 Cl 

"3 

z 

0 

N 

C3 

}} 

H 

A4  B4_ 

I 

y4 

A4  B7 

h 

0 I 

5ii 

E 


As  Bs  Ds] 


ys 

0 

L_ 


[' 


A,  Bc  D.  I 

5 5 


It  can  be  seen  that  from  point  5 the  first  covariances  between  the  coor- 
dinates are  involved. 


point  6 


S]  = V6+B6CS+D6C4  B6  D3 


I”  [\  »6 

51 


LC4J 


'CS 

. Fs  Bs  Df] 

~ys 

^ Bs  D- 

fo 

o 

1 

1 

C4 

0 I 0 

C3 
L U 

and  so  on. 

According  to  the  above  equations,  a new  gradiometer-aided  navigation  sys- 
tem space  traverse  is  carried  out.  The  results  are  listed  below. 
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— m m rvp 


_ AiT&r.m  TRACTly^'^ 


Point 

var.-cov.:  (m-) 

3 

0.1638  10"4  0.1438  10-4  0.2182  10'3 

0.1658  10'4  0.2182  10'3 
0.3314  10'2 

4 

0.8038  10'4  0.7190  ’O'4  0.1091  10'2 
0.8038  10'4  0.1091  1G"2 

0. 16S7  10'1 

5 

_T,  .3  -2  I 

0.3S43  10  • 0.3164  10  0.45SD  10  j 

0.3543  10‘3  0.4501  1C"2  ! 

0.7345  10*'  ! 

i 

6 

0.1514  10*  2 0.1352  10" 2 2'. 1971  10* 1 

-2  - ■ i 

0.15H  10  0.7051  10  1 

0.313' 

; 7 

I 

i 

0.6420  10'2  O.S'37  1C*2  0.8364  10' 1 i 

1 

0.642i'  10'2  0.S708  10' 1 ! 

| 

1.3315 

i 

From  the  given  result®  it  is  evident  that  the  omitted  covariances  do  not 
play  any  critical  role  on  the  derived  coordinate  variances. 

c)  Investigation  on  a multipoint  approximation: thi  last  hope  to  improve 
the  bad  behaviour  of  the  gradiometer-aided  inertial  navigation  system 
is  to  include  more  terms  for  the  Stirling’s  approximation  formula  and 
thus  to  approximate  better  the  inertial  acceleration  components . Eq. (5.2.5) 
includes  only  the  first  term  of  the  mentioned  formula,  but  manipulation 
of  the  second  term  gives  to  the  fundamental  .^nation  cf  inertiaj  naviga- 
tion to  be  simulated  the  following  form: 

Xi  ‘2At^AX.  , +GX  +GXX  (Xi-i  “Xi-2  )+gxy  (Vi  >♦ 

Ai-2,i-1  *i-2  AXi-2,i-l  11  12  ^i-2,1-1  11  1 ‘ 


♦G 


XZ, 


i_2>i_1  (Zi-1  -Zi-2  »*4Xi-2  *Xi-4 


(7.3.8) 
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y‘  ‘\-2.i-A-2  -V;.!.,  (X*-’  "1-1  -*.-a  >* 

<Z‘-'  ‘Zi-2  ),*4Vi-2  -"i-S  *Vi-> 

Z‘  (x‘-’  -X2-2  «t-i  -'.-2  >* 

*Szi-2.i-1  <Zi-’  _Z‘-2  ”*4Z‘-2  ',Zi-3  *Zi-4 


The  same  procedure  followed  by  the  previous  simulation  studies  is  applied 
again  to  the  above  equations.  The  space  traverse  for  the  multipoint  appro- 
ximation gives  the  following  results: 


From  the  above  listed  results  it  is  evident  that: 

1) the  navigation  system’s  behaviour  is  now  better 

2) the  tremendous  instability  of  the  Z-channei  has  been  already  diminishea. 

3)  the  system  gathers  less  errors  relative  to  the  previous  analysed  tra- 
verses but  it  still  needs  to  be  updated  at  a certain  navigation  time  after 
the  initial  observation  point.  This  comes  in  confrontment  with  all  up-to- 
date  aided  navigation  systems  which  should  be  filtered  out  continuously 

as  far  as  their  conrritted  errors  are  concerned.  In  the  majority  of  inertial 
navigation  applications,  Kalman  filtering  is  continuously  applied (say, eve- 
ry some  seconds  of  navigation  time)  to  reduce  the  systems’  inclination  to 
gather  errors  which  sometimes  are  intolerable. 


I 
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8.  Conclusions  and  Recommendations 


We  live  undoubtedly  in  the  space  age.  Everything  moves  quickly  and  accura- 
tely. Inertial  navigation  systems  have  been  already  introduced  to  accomp- 
lish the  latter  but,  for  the  time  being,  in  a limited  sense.  Well-known 
that  systems  which  can  guide  a moving  vehicle  are  burdened  with  errors (so- 
metimes very  big  ones)  and  it  is  depending  on  the  objective  of  the  mission 
that  the  guiding  system  can  be  considered  to  be  successful  or  not.  For  exam- 
ple, the  strict  requirements  imposed  on  a satellite  guidance  system  are  not 
applicable  in  a slow  moving  vehicle  such  as  a submarine.  Thus,  the  system’s 
efficiency  will  be  depending  on  the  objective  of  the  guided  object.  Under 
that  prism  we  a.ial/se  an  aided  navigation  system  in  case  of  terrestrial  nc  • 
V- gat ion  and  especially  when  the  moving  vehicle  is  a cruise  aircraft.  The 
external  aid  is  cons: seed  of  three  mutually  perpendicular  gradiometers,  a 
revoluntary  equipment  with  the  capability  to  measure  the  gravity  gradient 
field  of  the  space  in  which  the  navigation  takes  place.  Several  aids  car. 
be  introduced  instead  of  the  gradiometers  but  we  have  chosen  their,  in  order 
to  disperse  possiole  fears  about  their  usefulness  of  operation  in  the  inve- 
stigated application. 


The  inertially  referenced  acceleration  of  a moving  vehicle  can  be  obtained 
by  adding  the  gravitational  and  non-grav: rational  (or  apparent)3Cceleration 
and  neglecting  certain  small  terms  which  have  been  already  discusses  pre- 
viously. Approximating  the  second  inertial  derivative  with  the  Stirling’s 
formula  and  then  solving  the  resulting  equations  with  respect  to  the  un- 
known coordinates,  the  error  propagation  law  can  be  applied  through  in 
order  to  get  the  expressions  to  be  simulated.  The  analysis  proves  that  a 
gradiometer-aided  inertial  navigation  system  is  very  unstable  like  all  na- 
vigation systems.  Particularly,  the  Z-channel  gathers  the  biggest  ammount 
of  errors  relative  to  the  rest  two  ones  and  after  some  seconds  th*  naviga- 
tion system,  as  an  instrumental  package,  collapses  as  far  as  its  perfor- 
mance is  concerned.  Aiding  the  system  with  a bare... o'er  or  an  altimectv 
or  generally  with  an  instrument  which  can  produce  any  kind  ol  height  iiur'or- 
mation  of  the  instantaneous  position  of  the  moving  vehicle,  then  the  insta- 
bility of  the  discussed  channel  can  be  effectively  reduced.  Relevant  to  the 
two  parameters  of  great  interest,  namely  the  acceleration  and  gravity  ones, 
the  former  causes  to  the  system  an  error  up  to  the  order  of  7% (considering 
its  contribution  for  some  first  points)  and  the  latter  a very  small  percen- 
tage of  error.  These  results  can  justify  the  conclusion  that  gravity  gradi- 
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ometers  perform  or  rather  behave  excellent  on-board  a moving  vehicle,  such 
as  an  aircraft,  or  that  the  gravity-induced  position  error  is  a negligible 
small  quantity. 

Considering  detailed  error  models  for  acceleration  and  gravity  gradient  mea- 
surements, a new  statistical  analysis  proves  that  the  aided  navigation 
system  under  investigation  remains  unstable  in  almost  the  same  fashion,  acce- 
lerometer measurements  errors  are  reduced  down  to  101  with  respect  to  the 
previously  derived  numbers  and  gradiometers  continue  to  fit  excellent  on- 
board. Trying  to  go  out  of  this  undesired  "cul-de-sac"  three  new  investi- 
gations are  carried  out: 

a) Making  the  operation  of  the  system  as  fine  as  possible,  then  the  behaviour 
of  the  navigation  system  does  not  change  appreciably.  It  will  break  down 
sooner  or  later. 

b) In  both  simulation  studies  the  initial  value  problem  of  inertial  naviga- 
tion is  under  consideration.  The  big  influence  of  the  two  initial  sets  of 
coordinates  on  the  derived  accuracy  of  the  system’s  output,  justifies  the 
hope  that  something  could  lie  inside  the  omitted  covariances  between  the 
coordinates  of  the  two  first  points.  The  statistical  analysis  shows  that 
the  navigation  system  is  again  dying  out  in  exactly  the  same  manner  as 
previously. 

c)  Since  the  whole  analysis  is  based  upon  the  approximation  made,  namely 
the  Stirling’s  formula,  it  is  reasonable  to  investigate  the  case  in  which 
more  terms  are  taken  into  account.  The  simulated  equations  change  format 
including  now  more  terms (or  more  initial  points).  The  new  statistical  ana- 
lysis justifies  the  hopes  that  the  system  cannot  be  only  badly-behaved  but 
it  possesses  the  ability  to  be  more  accurate  or  for  that  matter  usable. 

The  Z-channel  instability  is  decreased  effectively,  the  system  gathers 
less  errors  than  previously  and  it  can  now  be  seen  with  hope. 

In  spite  of  the  not  so  accurate  approximation  for  the  second  inertial  de- 
rivatives, the  quantization  error  to  be  committed  is  a very  small  quantity 
relative  to  the  total  error  budget.  Trying  to  give  an  order  of  magnitude 


for  that  error, I could  simply  mention  that  if  the  navigation  system  ope- 
rates every  O.lsec.,  then  after  20  minutes  of  flight  the  quantization  error 
is  up  to  the  ordeT  of  4cm. 

I i 

Further  studies  on  that  subject  are  advisable  according  to  the  following 

L 


general  ideas: 

1) Try  to  include  more  terms  in  approximating  the  second  inertial  derivati- 
ves in  deriving  the  simulated  equations  for  the  first  two  or  three  unknown 
points.  Simulate  them  and  revise  them  for  the  next  few  points  but  include 
now  more  terms  as  far  as  the  Stirling's  formula  is  concerned.  Apply  this 
procedure  continuously  until  the  system  has  gathered  errors  which  cannot 

be  further  tolerated.  Then  a filtering  technique  can"refresh”  the  systar. 
in  order  to  begin  again  the  above  discussed  procedure.  In  view  of  the  obtai 
ned  results  of  the  multipoint  approximation  studies,  I strongly  believe 
that  the  proposed  analysis  will  turn  out  to  be  very  fruitful. 

2) Apply  the  well-known  technique  of  Kalman  filtering.  Extensive  literature 
addresses  this  problem  and  therefore  it  will  not  be  discussed  herein. 

Perhaps  one  might  oe  keen  on  asking  why  I insist  so  nuch  on  applying  gravi- 
ty gradiometer  techniques  for  terrestrial  navigation  systems  since  the  up- 
to-date  used  systems  can  operate  with  better  accuracy (one  nautical  mile  per 
hour  flight,  approximately).  The  detail  and  accurate  notion  of  the  earth’s 
gravity  field  is  not  of  great  importance  for  such  applications.  But  the 
answer  being  straightforward,  comes  with  the  question:  what  these  gradiometer 
unaided  systems  can  do  for  space  missions  in  cases  of  which  the  or.-uoard 
platform  travels  through  different,  successive  and  completely  unknown  gra- 
vity fields?  I believe,  they  can  do  nothing (strictly  speaking). 

A gradiometer-aided  inertial  navigation  system  turns  out  to  be  very  profi- 
table in  many  applications  and  it  can  surely  upgrade  the  abilities  of 
Mankind  to  explore  the  mystery  which  is  here  to  stay,  the  distant  ccsmcs. 
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Appendix  A:Gon?uter  Pro  grannies  Used 


Similation  I uni  II 


PROGRAM  DOOKAK 

DIMENSION  A < 1 9 f 3 1 r£><l9?19)  f C ( 3 ’ 3 ) f ZW  < 1 9 » 1) 

REAL *8  AfPfC 

rt--3 
MM--  19 
NN-19 

CALL  ALESEN(AfMfNfO> 

CALL  ALESEN' R.nh.NN.2) 

CALL  MAMUL3  < A f H • N f 1 > B f MM  f NN  tO?A<MrNtOiCtNiNi IW  f 1 > f 1 J 

CALL  AORUCK  (C.NfN.O) 

STOP 

END 

SUBROUTINE  MAMUL3 (AfHAfNAfKENNAf BfMB fNBfKENNBfC fMC fNC fKENNCf 
l R . MR . NR f ZU  f MZW  f GYM ) 

01  ME  NS  ION  A < MA  • NA .MMBfNB)  fC<mCfNC>  fR(MRfNR)  fZW<MZWf1> 
c’EAl  *8  A.tuO.r  AF  l • F AF  2 * F AN  3 f R • 3 * ZW 
INTEGER  AZ.ASfCZfLSfSYM 
LOGICAL  NENNAL*NENNBLfKENNCLfSYHL 
kennal =kenna . ne . i 

KENN»L*NENNB.NE. t 
•\ENNCL*KENNC.NE . I 
GYML-SYM.EO. 1 
AZ*-MA 

as-na 

IF  IKENNAL)  GOTO  41 
AZ-NA 
AS*MA 
•»l  CZ*MC 
CS*NC 

IF'NENNCL)  GOTO  42 
CZ *NC 

>:s  -mc 

♦2  HOI  NMfAZ 

002  I=1fCZ 

c=-o. 

ti02  J-1fAS 

IF(.NOT.NENNALi  FAK1=A(JfK) 

IF  iKENNAl  ) FAM»A<KfJ> 

IF (.NOT. NENNBL)  FAN2“B(IfJ> 

Ir (hENNBL ) FAK2»H<JfI) 

5*S»FAK1*FAK2 
2 ZU<I f 1 ) -S 
>1 

IF(SYML)  J-K 
001  L“JfCS 
S*0. 

001  I-1fCZ 

IF  < . NOT . KENNCL)  FAN3«C(LfI> 

IF  vKENNCL)  FAK3-C(IfL) 

'5«=S+ZW  (I  f 1 ) *FAK3 
R(KfL  >=S 

t IF'SYML)  R<LfK)f=S 
RETURN 
END 

SUBROUT I NE  AINVER  < A f N) 

SUBROUTINE  NATRITZENINVERTZERUNG 
DIMENSION  AvNfN) 

IF<N.NE.1>  00 TO  11 

A«l'l>»l/A(ltl> 

GOTO  12 
11  0010  K-2rN 
Kl-K-1 
0010  I-1fK1 
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XO  A'Kf l **0 
Nl*-N-  1 
ClOl  K*>lfNl 
T=A'KfK> 

A<KfI>>*X. 

K "-K  ♦ L 
602  I=K?»N 
?-«A  CK  r I > /T 
60?  J 1 r*  1 1 K 

1 A(If J1)=*A<IfJ1>-S«A<Kf Jl> 

604  J2«IfN 

A A<If J2)*A<If J2)-S*A(KfJ2> 

2 CONTINUE 
1 A(K.K)*T 

C INVERSE  MATRIX 

r«A  <ntN> 

604'  fv  ^ ' • N 
* A tNif.)  eft  (N  f K)  /T 

por.  xi^-ifni 
r«N-i j 
7=A(I .1) 

A'If I>*  X . 

606  K»1»I 
S=A* I *K> 

12* 1*1 

607  J*I2.N 

7 S*S-A(I. J)«A(J.K> 

/.  A ( 1 1 K)  *S/T 
UMSREXCHERN 
HOC  K=I2fN 
■■j  A'IfR)=A<RfI> 

12  RETURN 
FN6 

SUDRGUTI NE  AB6SUB  iNENN » A*  MA  » NA  f RENNA . 6 . M6 . Nt  * NE NNB f C f MC . NC  i 
C SUBROUTINE  MATRIT7ENA60ITI0N  UN6  -SUBTRAM  ION 

El  I ME  JSXON  A ( MA  • NA ) fB'HBfNB)  fC(MCfNC) 

INTEGER  A7*  AS 

LOGICAL  KENNLf KENNAL fKENNBL 

KENNL«KENN.EQ. 1 

KENNAl*KENNA.NE. 1 

KENNPL  *KENNB . NE . 1 

AZ^MA 

AS*NA 

IE  (KENNAL)  GOTO  41 

A7-NA 

AP*MA 

*1  IiOl  K*1fAZ 

601  I- If  AC 

IF  < . NOT . KENNAl  > FAK1*A(IfK> 

IF  < KENDAL)  FAK1»A<KfI) 

IF < .NOT.KENNBL)  FAK7“B<IfK> 
tF(KENNBL)  FAK?*BfKfI> 

C <K  f I i "FAUX  *FAK2 
1 IF.KENNL)  C (Kf 1 ) “FAK1-FAK2 
RETURN 
EN6 

SUBROUTINE  ABRUCK (AfMfNfKENN) 

C SUBROUTINE  BRUCKEN  EINER  MATRIX 
DIMENSION  A<«# N) 

EEAL*6  A 
Ni*N 
M1«M 

IF (KENN. NE. 1 ) GOTO  43 
H1*N 
N1*M 

*3  WRITE <6f101 -MIfNI 
101  FORMAT (215/ ) 

DO  1 I-1fM1 

IF (KENN. NE . 1 ) WRITE <Af 102)  (A • I f J) t J* 1 .NX > 

IF (KENN.EO. 1 ) WRITE  <A* 102)  (A i Jf I) f J*1 1 N1 ) 


•PBIlC* 


-aSSSS**' 


i** 


rj  (.« 
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tor  rORMATUH  * 3F  14.4) 
l WRITE (61 102) 

RETURN 

ENIi 

SUBROUTINE.  A L E b'  E N (A.M.N.KENN) 

•'  LESFN  FINER  UNTEREN  IiREIECKShATRIX 

DIMENSION  A<M,N> 

REAL«e  fli < 

M1*M 

N1=N 

IF(KENN.NE.l)  GOTO  SI 

M1*N 

N1»M 

SI  IF  <KENN.GE . 2) GOTO  32 

DO  1 1*1. Ml 

IF (KENN.EQ. 1 > READ <5. 102)  (A ( J, I > , JM ,N1 > 
IF  <KENN. EQ.  0)  READ  (3. 102)  (A  ( I . J)  . J=  1 » N1 ) 
102  FORMAT <8F 10. A) 

1 CONTINUE 
3 READ  <3. 102>  X 

RETURN 
DO  2 1*1.  M 

READ (S. 102) <A ( I . J) . J* 1 . I) 

IF (KENN .EG . 3> GOTO  3 
DO  A 1*1, M 
DO  A J* 1,1 
•4  A ( J , I ) -A  ( I , J> 

GOTO  3 
END 
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Programmes  for  the  Quantization  Error  Studies 


a. Variance -case 


PROGRAM  TWINP1 
REAL*€  A . b » C 

REAPU»50>  NN.tiTrSA.TA.SGf  TG 
SA2«SA*« J 

SG2-Sr.*»2  . 

CALL  UAF (PT.SA2.SG2. TA.TG.NN) 

SO  FORMAT (15. 3F 10. «. IF 15. 10. IF  10. a) 

STOE 

ENP 

5UPE0UT I NE  VAR  ( IiT  . S A? . SS2 . T A . TG . NN ) 

PO  5 IN*-3«NN 

X~0 

tiO  6 IM--2.INM 
i >0  c<  IR— r.lNh 

y=X*HFLOAT  t (IN-Ih'  « < IN-  IR)  ) »PEXE  (-PAPS  <PT*l'f  Lu'  TM-  IR)  ) /TA) 
o CONTINUE 
Y-O 

PO  7 IM-2.TNN 
t'O  7 IR-.'«  INM 
PO  7 INi.;?,  in 
PO  7 IE--O.IR 

Y^r-t  riFtOAT  ( »IN-In)  « (IN-IR)  ) « PEXE  (-PAPS  (P I *PK  LOA  f ' IK-IF  ' ) / TG) 
7 CONTINUE 

S=PT**4#  < ?.A2«X  + SC 2* Y ) 

WRITE (6.1 00)  IN-? 

WRITE (2.100)  IN.? 

! TO  FOEnAT  <10X»F1j.m> 

RETURN 

fNTi 


b. Covariance  case 


■•■•ROGFAn  TUINS3 

REAL  *8  Cir.SAf  TA.SG»TG.SA2f3G?.S 
DIMENSION  S(20f.:0' 

F'EMD  1 ] * 50  1 NN  f D T • CA  f TA  f 3l3  f TG 
EA^SA*# ,? 

CCS  3G**2 

CALL  COW  < fiT  f SA2  • 3G2  f T A f TG  f NN . S> 

NNM?»NN~2 

50  FOR'hA T < 1 5 f 3F10.4. 1F15.  10f  1F10.4) 

urite  ( i?,f  :oo* 

WRITE  \e>*20O) 

WRITF  2«200) 

200  FORMAT  (lHt  f 12H  VAR  - COV  =) 

DO  a I *•  1 f NNM2 

JRITC  <13f 100)  \3iIfJ>  fJ»1fNNMC> 

JRTTE  . : F j 00 F I F .1)  . J*.  1 F NNM2) 

■fRJte  ioo>  is-  r*  j>  .j-Lfnnmzo 

100  FORMAT  < 5>f12(D10.2)  ) 

CONTINUc. 

STOF 

END 

S'  TROUT I NE  COVA  - riT  . SAC  f SG2  f TA  f TG  f NN . S > 

DIMENSION  3<20f 20  > 

REAL  *8  DT.SA. TAfSGf TGfSA2fSG2f Sf X f T 
DO  5 IN-3  f NN 
DO  5 IS^SfNN 
*■0. 

INM- IN-J 
ISM-IS- 1 
DO  6 TMCfINM 
DO  a IR*  _>»  ISM 

X=X+DFL0AT . (IN-th>* (IS- IR) > ffDEXF  (-DADS (DT«DFLOAT \IM-IR) ) /TA) 

4 CONTINUE 

Vv,>. 

DO  7 I M*-  2*1  NM 
DO  ’ TR-’-ISM 
DO  7 IT  =2 f IM 
DO  7 IF»2fIR 

f * Y+DFLOAT  < < IN  -IM)  • ( IS- IR)  > «DEXF'  (-DADS  (DT*l'F  LUA  T ‘ IK- IF')  ' /TG' 
7 CONTINUE 

S I IN-2f IS -2) »DT- •-* iSA2*X+SG2»Y) 

5 CONTINUE 
RETURN 
END 


